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PREFACE 


This  paper  is  dedicated  to  the  memory  of  my  son,  Brian  Jon  Elchlnger. 

The  purpose  of  this  study  was  to  find  a  method  of  predicting  the 
environment  In  the  path  of  a  Mach  stem  In  a  nuclear  explosion.  There  are 
many  Indications  of  the  Inadequacy  of  present  modeling  methods.  The 
prediction  method  and  program  In  Appendix  F  Is  Intended  to  be  simple,  but 
accurate.  As  the  study  progressed  It  became  obvious  that  the  Mach  stem 
phenomenon  Is  a  dynamic  process  that  changes  throughout  Its  life.  Thus 
to  adequately  describe  the  conditions  would  require  an  Inordinate  amount 
of'' equations.  This  study  Is  self-limited  to  the  region  of  Interest  of 
survivability  studies  of  aircraft. 

In  performing  this  study,  I  am  deeply  Indebted  to  Major  Larry  McKee, 
my  advisor,  for  his  assistance  and  guidance.  I  also  wish  to  thank 
Dr.  George  Ulrich  of  the  Defense  Nuclear  Agency  and  Mr.  Ray  Ruetnlck  for 
their  advice  and  help.  Finally,  I  wish  to  thank  my  wife  Leanne  and 
daughter  Heidi  for  their  patience  and  understanding  during  this  period. 


William  E.  Elchlnger 
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Abstract 


'  ^  k  semlemplrlcal  model  was  developed  for  the  treatment  of  the  Mach  stem 
region  of  a  nuclear  alrburst.  This  model  predicts  the  conditions  that 
would  be  observed  by  an  aircraft  or  missile  as  well  as  the  limits  of  this 
region  in  space.  The  model  is  based  upon  current  shock  theory.  Where 
three  shock  theory  falls  to  accurately  predict  physical  reality,  new 
relations  are  developed  or  empirical  data  is  used. 

Specifically,  this  model  predicts  the  path  of  the  triple  point,  the 
overpressure,  dynamic  pressure,  time  of  arrival,  and  direction  of  the 
shock  impulse  above  the  ground.  An  explanation  of  the  development  of 
each  prediction  is  made  and  compared  to  actual  nuclear  or  high  explosive 
test  data.  Additionally,  a  comparison  is  made  between  conventional  Mach 
stem  modeling  and  this  model. 

Unique  to  this  model  is  a  method  of  predicting  the  variation  of 

pressure  with  altitude  above  the  ground.  For  low  scaled  heights  of 

<r 

burst,  the  overpressure  found  at  an  altitude  of  20^percont  of  the  triple 
point  height  is  greater  than  that  on  the  ground.  In  addition,  the 
overpressure  measured  just  below  the  triple  point  is  found  to  be  only  60  ^ 
percent  of  the  ground  overpressure  scaled  to  altitude.  This  prediction 
is  radically  different  than  conventional  ground  overpressures  scaled  to 
atmospheric  pressure  at  altitude.  The  predictions  made  by  this  model  are 


verified  by  nuclear  and  high  explosive  test  data. 


NOTATION 


Air  Zero 
Ground  Zero 
HOB 
HTP 

Incident  Shock 
M 

OTP 

P 

dP 

Reflected  Shock 
Ta 

Triple  Point 

U 

a 

y 

6 

0 

p 

T 

0 

0 

Subscripts 
1,2, 3,4 


Point  of  detonation  of  the  weapon 
Ground  location  directly  below  the  burst 
Height  of  the  burst  above  the  ground 
Height  of  the  triple  point  above  the  ground 
Shock  wave  generated  directly  from  a  burst 
Mach  number 

Origin  of  the  triple  point,  the  ground  distance  at 
which  the  Mach  shock  begins  to  form 

Pressure  (psl.) 

Overpressure 

Shock  wave  generated  by  reflection  of  the  Incident 

Time  of  arrival  of  the  shock  wave  (sec.) 

Point  of  Intersection  of  the  Incident,  Mach  and 
reflected  shocks. 

Velocity 

Angle  of  Incidence 

Ratio  of  specific  heats — equal  to  1.4  for  air 
Angle  of  deflection  of  gas  through  a  shock 
Hedge  angle 

Shock  strength  as  defined  by  Po/(Po+  dP) 

Density 

Direction  of  the  gust  loading  of  airfoils 
Angle  between  flow  of  air  and  shock 
Triple  point  trajectory  angle 

Denotes  region  of  measurement 


I.  Introduction 


It  can  be  shown  that  air  bursts  are  most  effective  in  attacking  soft 
or  moderately  hard  targets  with  nuclear  weapons.  Further,  surface  bursts 
loft  great  quantities  of  dust  and  debris  into  the  atmosphere.  This 
material  can  be  exceedingly  damaging  to  subsequent  warheads  targeted  in 
the  vicinity  of  the  first  blast.  Thus  it  is  likely  that  airbursts  will 
be  used  against  many  military  and  civilian  targets.  This  paper  deals 
with  air  blast  phenomenon  in  the  Mach  stem  region  of  the  blast  wave 
emanating  from  a  nuclear  weapon.  The  Mach  stem  is  an  effect  that  occurs 
in  near  surface  bursts  (less  than  1000  'j  aled  Feet)  and  which  can  reach 
up  to  five  miles  in  height.  It  is  thus  important  in  survivability 
studies  of  airborne  vehicles  such  as  aircraft  and  launching  missiles. 

Presently,  the  overpressure  in  the  Mach  stem  is  assumed  to  be  constant 
vertically  and  is  scaled  only  by  the  variation  of  pressure  due  to  the 
altitude  of  the  target  (19:2-50).  It  is  further  assumed  that  the  Mach 
stem  is  vertical  to  the  ground  at  all  altitudes.  This  leads  to  the 
assumption  that  the  Impulse  which  is  delivered  to  a  target  in  the  Mach 
stem  is  parallel  to  the  ground.  The  time  at  which  the  Mach  stem  blast 
wave  arrives  at  a  target  above  the  ground  is  not  addressed  at  all.  A 
cursory  examination  of  photographs  of  near  surface  nuclear  bursts  clearly 
shows  that  the  Mach  stem  is  in  general  not  vertical.  Changes  in 
curvature  in  the  Mach  stem  are  an  indication  that  the  overpressure  in  the 
Mach  stem  is  not  constant  over  the  entire  vertical  distance  of  the  shock. 

There  are  several  works  that  treat  the  phenomena  of  nuclear  bursts. 


The  most  important  of  these  are  the  Defense  Nuclear  Agency  manual  EM-1 , 
Glasstone's  Effects  of  Nuclear  Weapons,  and  the  Defense  Atomic  Support 
Agency  manuals  DASA-1200  and  DASA-2506.  Amongst  these  manuals,  blast 
effects  are  treated  in  great  detail.  However,  the  conditions  that  are  to 
be  found  in  the  path  of  the  Mach  stem  above  the  surface  of  the  earth  are 
not  discussed.  In  order  to  develop  a  model  which  will  predict  blast 
parameters  this  paper  will: 

1.  Examine  current  theory  of  shock  reflection, 

2.  Use  the  results  of  current  theory  to  predict  the  path  of  the 
triple  point. 

3.  Examine  nuclear  blast  data  to  find  relations  for  the  shape, 
time  of  arrival,  and  overpressure  of  the  shock  at  altitude. 

4.  Compare  predictions  made  by  the  resulting  relations  to  other 
predictions  and  experimental  data. 
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II.  Mach  Effect 


The  Mach  stem  was  first  discovered  in  1878  by  the  pioneer  researcher 
Ernst  Mach  for  whom  it  was  named.  It  remained  obscure  until  1941  when 
Von  Neumann  introduced  the  first  analytical  criterion  on  hov/  the 
phenomenon  occurs.  During  the  Second  World  War,  American  scientists 
studied  the  effect  and  discovered  that  the  destructive  effects  of  bombs 
could  be  increased  by  exploding  the  weapon  at  some  height  above  the 
ground.  The  blast  overpressure  as  measured  at  a  fixed  location  on  the 
ground  increased  as  the  height  of  detonation  increased  until  an  ideal 
height  was  reached  and  then  decreased  with  higher  heights  of  burst.  This 
increase  in  ove'rplressure  is  caused  by  the  formation  of  a  Mach  stem  at 
certain  ranges  and  overpressures  (17:1).  When  the  early  nuclear  tests 
were  done  in  the  mid-1950's,  investigations  were  conducted  to  determine 
if  the  same  effects  occurred  in  nuclear  weapons  as  well.  The  tests 
showed  that,  with  proper  scaling,  nuclear  weapons  exhibit  the  same  blast 
characteristics  as  conventional  weapons. 

When  a  nuclear  weapon  is  detonated  at  some  height  above  the  surface  of 
the  earth,  a  spherical  blast  wave  is  emitted  from  the  device.  This  wave 
strikes  the  ground  at  an  angle  of  Incidence  that  grows  from  zero  to  90 
degrees  for  locations  away  from  ground  zero.  The  angle  of  Incidence  is 
defined  as  the  angle  between  the  tangent  to  the  shock  and  the  ground. 

A  reflected  shock  wave  is  formed  which  travels  behind  the  incident 
wave.  Because  the  air  behind  the  incident  wave  is  hotter  and  more  dense 
than  the  air  in  front  of  it,  the  reflected  wave  travels  faster  than  the 
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incident  wave.  Under  certain  conditions,  the  reflected  wave  may  catch  up 
to  the  incident  wave,  merge  with  it  and  thus  reinforce  the  incident  shock 
wave  (See  figure  1).  The  region  in  which  there  is  a  separate  and 
distinct  reflected  shock  is  called  the  regular  reflection  region.  The 
region  where  the  shock  wave  is  reinforced  by  the  reflected  wave  is  called 
the  irregular  or  Mach  reflection  region.  The  top  of  the  Mach  stem  is 
called  the  triple  point.  At  this  point,  three  shocks  (the  Mach  stem,  the 
incident  wave,  and  the  reflected  wave)  converge. 


Figure  1.  Mach  Stem  Formation 

Several  other  effects  have  been  observed  in  the  nuclear  weapons  tests. 
For  low  scaled  heights  of  burst  (less  than  200  feet),  a  substantial 
portion  of  the  reflected  wave  must  travel  through  the  region  of  the 
expanding  fireball.  The  air  in  this  area  has  been  heated  to  extremely 
hot  temperatures  (7000  degrees  K.)  (27:307).  This  causes  the  relected 
shock  wave  to  be  refracted  as  it  passes  through  the  region.  The  net 
effect  of  this  refraction  is  to  weaken  the  reflected  shock  and  to 


Increase  its  velocity  in  the  vertical  direction  (See  figure  2).  The 
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Figure  2.  Shock  Wave  Contours 
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Figure  3.  Precursor  Structure 


velocity  of  the  reflected  wave  may  be  Increased  to  the  point  that  the 
part  which  passed  through  the  fireball  overtakes  the  incident  wave. 

Under  normal  conditions,  the  reflected  wave  would  overtake  the  incident 
wave  close  to  the  ground  first,  forming  a  Mach  stem.  For  low  heights  of 
burst,  the  reflected  wave  above  the  fireball  overtakes  the  incident  wave 
and  forms  a  second  Mach  stem  directly  over  the  burst.  The  size  of  this 
region  of  coalescence  can  be  estimated  using  techniques  described  in 
DASA-1200  (22:4-177).  However,  since  the  reflected  wave  is  greatly 
weakened  by  the  passage  through  the  fireball,  the  effect  of  its 
reinforcement  of  the  incident  shock  is  minimal.  In  the  Teapot  series  of 
tests,  for  example,  this  closure  was  seen  on  two  of  the  tests,  numbers  4 
and  12.  On  shot  12,  there  was  no  perceptible  reinforcement  of  the 
incident  wave  at  all  in  the  region  directly  above  the  burst.  On  shot  4, 
the  coalescence  of  the  waves  corresponded  to  an  increase  in  effective 
yield  of  only  16  percent  (46:42-45,76-79).  For  an  ideal  Mach  stem  formed 
above  a  perfectly  reflecting  plane,  the  increase  in  effective  yield  would 
be  100  percent.  Because  the  effect  of  the  closure  of  the  shock  waves 
above  the  burst  does  not  greatly  increase  the  threat  to  aircraft  it  will 
not  be  examined  in  this  paper. 

The  radiation  emanating  from  a  nuclear  weapon  can  also  cause  a 
phenomenon  known  as  a  precursor.  This  radiation  heats  the  air  near  the 
ground,  causing  the  incident  wave  to  travel  faster.  A  toe-shaped 
protuberance  is  formed  on  the  incident  shock  wave  near  the  ground.  This 
is  called  the  precursor  (See  figure  3).  This  effect  is  very  complex  and 
not  well  understood.  There  are  great  variations  in  overpressure  with 
height  in  the  precursor.  This  makes  an  accurate  estimate  of  the  actual 
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overpressure  on  the  ground  very  difficult.  The  angle  of  incidence  of  the 
shock  wave  in  the  precursor  region  also  changes  dramatically  (55:96-97). 
The  precursor  effect  can  mask  the  formation  of  the  Mach  stem  in  its  early 
stages.  For  this  reason,  the  actual  origin  of  the  triple  point  (the 
point  at  which  the  Mach  stem  begins  to  form)  is  not  well  known  for  many 
heights  of  burst.  Because  of  the  uncertainty  of  the  conditions  in  the 
precursor  region,  it  will  not  be  considered  in  this  study. 

In  recent  years,  laboratory  experiments  have  discovered  that  there  are 
several  different  forms  of  Mach  reflection,  namely  single  Mach 
reflection,  complex  Mach  reflection,  and  double  Mach  reflection.  The 
phenomenon  observed  at  large  distances  (about  equal  to  the  HOB)  is  the 
single  Mach  reflection.  At  very  close  distances  and  high  overpressures, 
complex  Mach  reflection  and  double  Mach  reflection  may  occur.  These 
effects  further  Increase  the  overpressure  on  the  ground  and  thus  may 
Increase  the  danger  to  reinforced  structures  on  the  ground.  Because  this 
report  primarily  concerns  the  effects  seen  by  airborne  vehicles,  the 
single  Mach  reflection  will  be  assumed  in  all  cases. 


III.  Oblique  Shock  Reflection  Theory. 

Use  of  Planar  Shock  Theory  for  Spherical  Shocks. 

Solving  the  fluid  dynamic  equations  In  a  spherical  geometry  Is  an 
extremely  difficult  task  both  analytically  and  numerically.  In  order  to 
simplify  the  problem,  a  spherical  shock  system  Is  assumed  to  be 
Instantaneously  modeled  as  a  planar  system  with  a  similar  configuration. 
Thus  the  spherical  shock  system  In  figure  4  Is  assumed  to  be  e 

Instantaneously  equivalent  to  the  planar  system  In  figure  5.  The  angle 
of  Incidence  of  the  spherical  shock  as  measured  by  the  tangent  to  the 
shock  Is  the  same  as  the  angle  of  Incidence  In  the  planar  shock.  The 
wedge  angle  9w  of  the  planar  shock  system  Is  equivalent  to  the  angle  B 
In  figure  4.  Once  the  pressure  and  triple  point  trajectory  angle  are 
determined  for  the  planar  system,  they  are  assumed  to  be  the  same  for  the 
spherical  system. 

One  series  of  calculations  Is  sufficient  to  determine  the 
configuration  of  the  planar  system.  Because  the  configuration  of  the 
spherical  system  changes  with  time,  it  requires  that  many  planar  systems 
be  solved  In  order  to  determine  the  characteristics  of  the  spherical 
shocks  with  time.  Solving  these  multiple  problems  will  enable  the 
prediction  of  the  path  of  the  triple  point. 

In  order  to  formulate  and  solve  the  problem  of  oblique  shock 


reflection,  several  assumptions  were  made.  The  medium  In  which  the  shock 
flows  Is  a  perfect  diatomic  gas.  For  this  gas,  gamma,  the  ratio  of 
specific  heats  is  equal  to  1.4.  This  restriction  can  be  later  relaxed  by 
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Figure  4.  Spherical  Shock  Configuration 


assuming  that  gamma  is  a  variable  function  of  the  incident  overpressure. 
The  second  assumption  is  that  all  fluid  flows  are  two-dimensional  and 
inviscid.  Flows  are  also  assumed  to  be  self-similar.  That  is,  the  three 
variables  x,  y,  and  t  are  not  independent.  The  problem  can  be  described 
at  any  point  in  terms  of  the  new  variables  x/t  and  y/t.  Thus  the 
components  of  velocity  in  the  x  and  y  directions  are  constant.  It  is 
further  assumed  that  when  two  shock  solutions  are  possible,  the  weaker  of 
the  two  solutions  will  be  the  physically  correct  one.  Although  it  has 
never  been  proven,  this  experimental  fact  is  assumed  to  be  due  to  a 
requirement  of  minimum  entropy  (5:1,2). 


Figure  6.  Regular  Reflection 
As  Viewed  In  the  Frame  of  Point  P  Stationary 
(Incident  Shock  moving  to  the  right) 


Two  Shock  Theory. 


Two  shock  theory  deals  with  the  reflection  of  an  oblique  shock  wave 
from  an  Ideal  reflector.  When  an  oblique  shock  wave  strikes  a  surface,  a 


reflected  wave  is  formed  (See  figure  6).  The  strength  and  angle  of 
reflection  may  be  found  from  the  analysis  of  the  air  flow  in  the  vicinity 
of  the  point  of  reflection  (Designated  P  in  figure  6).  With  respect  to 
an  observer  traveling  on  the  point  of  reflection,  P,  the  ambient  air 
(region  1)  is  flowing  towards  him  with  a  velocity  equal  to  Us/sinOf.  The 
flow  of  air  is  depicted  in  figure  6  by  arrows.  The  conditions  across  the 
shock  from  region  1  to  region  2  are  expressed  by  the  following  four 
equations  which  are  statements  of  conservation  of  mass,  momentum  and 
energy:  (5:8,9) 


Pjtana=  p2taTiiOt  -  (III-l) 

PlU^sina=  P2U2sln(a-  (III-2) 

?!+  PlUi^sin^a  =P2+  /^U2*sin“(a  -  di)  (III-3) 

H^+ljUj^sin^a  =H2+i5U2*sin*(a  -  dj^)  (III-4) 


Where:  are  the  air  densities  in  regions  1  and  2. 

^1,U2  are  the  velocities  of  air  in  regions  1  and  2 

Us  is  the  velocity  of  the  incident  shock 

^1,^2  pressures  in  regions  1  and  2. 

Otis  the  angle  of  incidence. 

dj  is  the  angle  of  deflection  of  the  incoming  flow. 

Hj,H2  are  the  enthalpies  of  the  air  in  regions  1  and  2. 

A  similar  set  of  four  equations  can  be  written  for  the  conditions 
across  the  reflected  shock  in  going  from  region  2  to  region  3.  Because 
the  flow  near  the  wall  must  move  parallel  to  the  ground,  the  strength  and 
angle  of  the  reflected  shock  must  be  such  that  it  causes  a  deflection  of 
the  flow  through  an  angle  equal  In  magnitude  to  the  first  deflection  but 


opposite  in  direction.  That  is: 


62  "  -  h  (III-5) 

These  nine  equations  and  the  equation  of  state  are  sufficient  to 
completely  describe  the  conditions.  Knowing  the  incident  angle,  the 
ambient  air  conditions,  and  the  shock  strength,  one  can  solve  for  the 
reflected  shock  (5:9).  As  the  angle  of  incidence  increases,  the 
magnitude  of  the  reflected  shock  Increases  until  at  some  angle  the 
strength  of  the  reflected  shock  is  greater  than  that  of  the  incident 
shock.  Thus  for  an  aircraft  in  this  region,  the  reflected  shock  will  be 
more  damaging  than  the  incident  shock.  The  angle  at  which  the  reflected 
shock  becomes  greater  than  the  incident  shock  is  found  from  the  relation: 

a  =^arccos((  y-l)/2]  (III-6) 

For  ideal  air  with  gamma  equal  to  1.4,  this  angle  is  39.23  degrees 
(15:327).  As  the  angle  of  Incidence  increases  past  this  angle,  the 
reflected  shock  Increases  in  magnitude.  In  solving  the  above  equations, 
one  will  find  that  there  is  no  real  solution  to  the  problem  for  certain 
shock  strengths  and  angles.  These  are  the  regions  of  Mach  reflection. 

The  angle  at  which  Mach  reflection  begins  for  a  given  shock  strength  is 
called  the  extreme  angle  and  is  designated  Qe.  This  angle  can  be  found 
from  the  following  relation  (29:54): 

[2x(  »7*x*+l)]*  -  4[[  y(  17  (l+x*(  ?7-l ))-!)+  »?(l+x*)]*-[  I7*x*+1]*]  «  0 


Where  x”  tan  (5r/2  -  Oe) 


(III-7) 


y+i  +  y-i  j  +  y+ij 


The  origin  of  the  triple  point  (OTP)  is  related  to  0(e  by  the 
following: 


OTP  »  HOBtanOte 


(III-8) 


Figure  7  is  a  graph  of  Ote  versus  the  incident  shock  strength. 

Plotted  as  shown,  the  curve  describes  the  boundary  between  the  Mach 
reflection  region  and  the  regular  reflection  region.  It  should  be  noted 
that  experiments  have  consistently  shown  that  regular  reflection  persists 
several  degrees  beyond  that  which  is  predicted  by  this  theory.  The 
difference  in  predictions  is  less  than  four  degrees  except  for  shocks  of 
strength  <.2  (29:55).  Shock  strength  is  defined  as  the  ratio  of  the 
pressure  ahead  of  the  shock  to  that  behind  it.  Thus  strength  is  always  a 
number  less  than  1  and  the  stronger  the  shock,  the  lower  the  'strength'. 
Note  that  the  TNT  test  results  shown  on  figure  7  are  consistently  several 
degrees  higher  than  the  predicted  angles.  Experiments  with  planar  shock 
waves  on  wedges  show  results  similar  to  the  TNT  test  results  (6:592). 

The  persistence  of  regular  reflection  past  the  predicted  point  is  known 
as  the  Von  Neumann  paradox.  Many  theories  have  been  advanced  to  explain 
the  persistence,  but  none  have  gained  widespread  acceptance  (61:180), 
(6:600-2). 

In  contrast  to  the  TNT  test  results,  the  data  from  the  nuclear  tests 
does  not  show  any  particular  pattern  except  that  the  triple  point  forms 
before  the  theory  predicts.  The  precursor  is  believed  to  be  the  reason 
for  this.  The  radiation  from  the  weapon  heats  the  air  near  the  ground 
and  changes  the  properties  of  that  air.  Thus  the  Mach  stem  may  form 
earlier  and  rise  faster  in  this  region.  When  the  precursor  effect  dies 


out,  the  path  of  the  triple  point  becomes  very  irregular.  It  has  even 
been  observed  to  decrease  in  height  as  the  precursor  died  out  and  then 
increase  again  along  a  path  closer  to  predictions  based  upon  theory  and 
HE  experiments.  Because  of  the  amount  of  dust  near  the  ground  and  the 
precursor  effect,  the  actual  origins  of  the  triple  points  are  not 
accurately  known.  Much  of  the  test  data  gives  conflicting  results  and 
thus  many  conflicting  data  points  can  be  quoted  from  the  literature. 
Shock  photography  does  not  give  good  results  until  the  Mach  stem  has 
risen  above  the  dust  and  ground  clutter  (which  is  generally  a  height  of 
at  least  ten  feet).  For  a  scaled  height  of  burst  of  400  feet,  the  Mach 
stem  may  have  traveled  over  2500  feet  by  the  time  it  gets  to  these 
heights.  Thus  backwards  extrapolation  to  find  the  OTP  is  not  very 
reliable  (See  also  (45:195)). 


Three  Shock  Theory. 

Figure  8  is  a  diagram  of  the  notation  to  be  used  in  this  discussion  of 
three  shock  theory.  There  are  four  general  regions  to  be  considered  in 
the  analysis  of  the  three  shock  configuration.  Region  1  is  the  ambient 
air  and  region  2  is  the  shocked  area  between  the  incident  and  reflected 
waves.  Regions  3  and  4  are  the  area  behind  the  Mach  shock.  The  air  in 
region  3  has  passed  through  both  the  incident  and  reflected  shock  while 
that  in  region  4  has  passed  through  only  the  Mach  shock.  The 
configuration  will  be  analyzed  from  the  point  of  view  of  an  observer 
traveling  with  the  triple  point.  To  the  observer  in  the  laboratory  frame 
of  reference,  the  triple  point  rises  at  a  constant  rate  which  can  be 
expressed  as  the  angle  As  stated  before,  the  flows  are  self-similar. 

That  is,  the  dimensions  of  the  flow  are  a  function  of  time  only. 


Although  the  Mach  stem  increases  in  height  with  distance,  the  angle  ^  is 
constant  throughout  the  experiment. 

In  the  laboratory,  shock  tubes  and  wedges  are  used  to  create 
nonstationary  Mach  reflections.  The  incident  shock  is  created 
perpendicular  to  the  wall  of  the  tube.  A  wedge  of  some  angle  6w  is 
mounted  in  this  tube.  When  the  incident  shock  strikes  the  wedge,  it  has 
an  effective  angle  of  incidence  equal  to  90  -  ^w.  Using  wedges  of 
different  angles,  the  point  of  transition  to  Mach  reflection  can  be  found 
and  plots  of  the  angle  l/f  versus  wedge  angle  and  shock  strength  can  be 
made.  Figure  8  shows  the  system  as  it  would  be  seen  through  a  window  in 
the  tube. 


Reflected  shock 


Region  2 


Incident  Shock 


Region  3 
Slipstream 


Region  1 


Mach  Shock 


Region 


Wedge  angle 


Figure  8.  Mach  Stem  Configuration 


To  an  observer  on  the  triple  point,  the  ambient  air  is  flowing  towards 
him  at  a  speed  of  =  Us/sin  where  Us  is  the  incident  shock  velocity 

and  <j>^  angle  shown  in  figure  6.  0]^  is  equal  to  90-^-  The 

velocity  of  this  ambient  air  is  parallel  to  the  trace  of  the  triple 
point.  In  describing  the  three  shock  configuration,  twelve  equations  may 
be  written.  They  are  statements  of  conservation  of  mass,  momentum  and 
energy  across  each  of  the  three  shocks. 


Pj^tan 

= 

P2tan(  01-  ii) 

(III-9) 

Pj^Uj^sin 

= 

P2U2Sln(  01+ 

(III-IO) 

^1+  PiUj^sin" 

= 

^2"*"  ^1^ 

(III-ll) 

= 

H2+^U2=*sln^(  01-  ap 

(III-12) 

/l^tan  <t>2 

= 

P3tan(02-  62) 

(III-13) 

P2U2sin  ^2 

= 

P3U3Sin(  02“  62) 

(III-14) 

^2+  /^U2*sin*  <t>2 

= 

P3+  P3U3“sin*(  02-  62'> 

(III-15) 

H2+^U2*sin®  ^2 

= 

H3+^U3*sin*(  02-  ^2) 

(III-16) 

Pi  tan  03 

= 

P^tanC  03-  ^3) 

(III-17) 

PlUjsin  03 

= 

P4U4Sin(  03-  ^3) 

(III-18) 

?!+  PiUj^sin*  03 

P4+  P4U4*sin*(  03-  a3) 

(III-19) 

^l+^Ui*sin®  03 

= 

H4+^U4*sln*(  03-  ^3)  (5:14) 

(III-20) 

In  addition  to  these  general  relations,  several  boundary  conditions 
apply  that  are  related  to  the  existence  of  the  slipstream.  The 
slipstream  is  a  surface  that  separates  two  regions  of  equal  pressure  and 
direction  of  flow.  The  f.^o  regions  have  different  temperatures,  entropy 
and  density.  This  is  because  the  air  in  region  3  has  been  shocked  twice, 
by  both  the  incident  and  reflected  shocks,  while  the  air  in  region  4  has 
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been  shocked  only  once.  Because  of  the  equality  of  pressure  across  the 
slipstream,  the  boundary  condition  may  be  written: 

P4  =  P3  (III-21) 

The  direction  of  the  flow  along  the  slipstream  must  also  be  the  same 
on  both  sides.  This  implies  that  the  air  passing  through  the  Mach  shock 
must  be  deflected  through  the  same  angle  as  the  air  passing  through  the 
other  two  shocks.  Thus: 

^3  =  dj  +  ^2  (III-22) 

With  the  equation  of  state  to  relate  the  properties  of  enthalpy, 
temperature,  pressure  and  density,  the  above  14  equations  contain  IS 
variables  Uj,  Uj,  U3,  ,  P2.  P3,  P4,  T^,  T2,  T3,  T4, 

il,  62,  and  ^3  (5:15)).  In  setting  up  an  experiment,  the  only  Initial 
conditions  are  Pj^^  wedge  angle  is  also  known,  but  its 

use  to  define  any  of  the  other  angles  Involves  another  variable,  ^  .  In 
contrast  to  Two-Shock  theory,  Three-Shock  theory  is  indeterminate.  In 
order  to  solve  these  equations,  an  additional  assumption  must  be  made. 

The  most  common  is  to  set  the  angle  ^  equal  to  90—4^,  This  assumption 
is  equivalent  to  assuming  that  the  Mach  stem  is  straight  and 
perpendicular  to  the  wedge.  With  this  additional  relation,  the  system 
may  be  solved,  usually  by  some  Iterative  method. 

The  above  relations  describe  the  shape  and  motion  of  the  three-shock 
system  only  for  the  case  of  strong  shocks.  Figure  9  compares  the 
predictions  of  two  and  three  shock  theory  for  the  case  of  a  weak  and  a 
strong  shock.  It  should  be  noted  that  two  shock  theory  is  good  for  weak 


shocks.  The  failure  of  three  shock  theory  is  often  attributed  to  its 
failure  to  account  for  viscous  and  heat  conduction  effects  (29:56), 


(29:57) 


IV.  Prediction  of  the  Path  of  the  Triple  Point. 


Prediction  by  Theory 

In  order  to  predict  the  path  of  the  triple  point  for  spherical  shocks, 
one  must  be  able  to  predict  the  value  of  is  the  instantaneous 

trajectory  angle  of  the  triple  point.  In  planar-  shocks,  this  angle  is 
constant.  In  the  case  of  spherical  shocks,  as  the  triple  point  rises, 
the  angle  of  incidence  changes  continuously,  and  thus  the  angle  at  which 
the  triple  point  is  instantaneously  rising  changes  accordingly.  The  path 
of  the  triple  point  is  given  by  the  following  relation; 


Height  of  Triple  Point  =*  /  tan  *1^  dr 


(lV-1) 


Thus  the  ability  to  trace  the  path  of  the  triple  poii\t  is  based  upon 
the  ability  to  accurately  predict  the  angle  ^  for  all  cases  under 
consideration. 

In  order  to  predict  this  angle,  three  different  cases  had  to  be 
investigated.  Planar  shock  theory  as  outlined  in  the  previous  section 
will  give  good  solutions  only  for  the  case  of  strong  shocks. 

Experimental  data  is  used  when  the  incident  shock  Mach  number  falls  below 
2.  However,  when  the  triple  point  rises  to  a  height  greater  than  the 
height  of  burst,  there  Is  no  theory  to  predict  the  results  nor  Is  there 
any  experimental  data.  This  Is  caused  by  configuration  of  the  shocks  In 
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which  the  incident  shock  does  not  precede  the  system.  In  this  region  the 
author  developed  a  model  to  predict  the  value  of 


Strong  Shock  Solution 

The  solution  for  the  strong  shock  case  is  based  upon  a  method  advanced 
by  Mirels  (44:bl-3).  The  original  method  was  worked  out  for  the  case  of 
infinitely  strong  shocks  only.  This  derivation  has  removed  this 
assumption  and  can  accurately  predict  the  angle  0  for  incident  Mach  wave 
shock  numbers  as  low  as  2. 

Six  basic  relations  are  required  to  solve  the  problem.  The  notation 
used  is  the  same  as  that  used  for  developing  the  three  shock  theory  and 
can  be  seen  in  figure  6.  The  first  of  these  relations  describes  an  angle 
X  such  that:  (44:b2) 

k=ct-/3-e  (iv-2) 

Once  again  the  oblique  shock  relations  are  used  to  find  the  ratios  of 
pressure  and  density  across  the  incident  shock.  The  frame  of  reference 
is  the  one  in  which  the  triple  point  is  stationary.  The  density  ratio  is 
denoted  as  K1  and  the  pressure  ratio  is  denoted  as  K2  for  later 
convenience  and  simplification  of  the  final  relations.  Thus  the  pressure 
and  density  ratios  can  be  found  from: 

P?  =  (  y+l)Mi^cos^(g  -f  =  K1  (IV-3) 

Pi  2+(  y-i)Mi*cos*(d  +{p) 

£2  *  2y  Mi*cos*(«  +  0)  -  K2  (IV-4) 

^  y+1  y+1 
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Figure  10.  Regions  and  angles  for  Strong  Shock  Relations, 


The  next  relation  is  the  ratio  of  the  airflow  in  region  2  that  is 

tangent  to  the  direction  of  motion  to  that  which  is  perpendicular  to  the 

direction  of  motion.  The  component  of  velocity  of  the  Incident  flow  that 
is  parallel  to  the  shock  front  is  not  altered  by  the  passage  of  the  air 

through  the  front.  Thus  the  tangential  velocity  in  region  1  is  the  same 

as  the  tangential  velocity  in  region  2.  The  velocity  component  that  is 
perpendicular  to  the  shock  front  is  decreased  in  magnitude  by  an  amount 
that  can  be  measured  by  the  density  ratio  which  has  been  designated  as 
Kl.  This  ratio  of  airflow  in  region  2  defines  the  angle  in  the 
following  manner. 


Vt2  =  cot  ^  Vnj  -  P,  8ln(  -  Kltan(tf  +  0)  (IV-5) 

^^^2  Vn2  Vnj  Pj  cosC  %  + 

Similarly  in  region  4  the  angle  d  defines  the  angle  of  the  flow  as 
does  in  region  2.  d is  also  the  angle  between  the  slipstream  and  the 
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Mach  shock.  In  front  of  this  region,  the  air  is  moving  into  the  Mach 
stem  at  an  angle  of  0.  Thus  the  ratios  of  densities  and  pressures  are 
defined  by  the  following  relations  and  designated  K3  and  K4. 

Pi.  =  (  Y+l)Vl^^cos^(  lif)  =  K3  fIV-6) 

Pi 

^4  =  2yMi"cos=^(  0)  -  2:zi  =  K4  (IV-7) 

■^1  (y+l)cos"(d y+1 

Notice  that  this  formula  assumes  that  the  Mach  stem  travels  at  the 
same  velocity  as  the  incident  shock.  To  an  observer  on  the  triple  point 
who  is  assumed  to  be  stationary,  the  ambient  air  appears  to  be  moving 
towards  him  with  constant  speed  and  direction.  This  velocity  is  defined 
by  the  velocity  of  the  incident  shock.  Since  the  shock  configuration  is 
assumed  to  be  self-similar,  the  air  must  approach  the  Mach  stem  at  the 
same  speed,  but  with  a  different  direction  than  the  air  approaching  the 
incident  shock.  By  assuming  that  the  angle  between  the  Mach  shock  and 
the  incoming  air  is  0,  one  is  assuming  that  the  Mach  shock  is 
perpendicular  to  the  wall  and  remains  that  way. 


Figure  11.  Shock  Configuration  in  Frame  of  Triple  Point  Stationary 


Figure  II  Is  a  simple  diagram  of  the  shock  configuration  from  the 
point  of  view  of  an  observer  riding  on  the  triple  point.  This  observer 
sees  three  stationary  shocks.  This  observer  also  sees  that  the 
overpressure  In  region  4  Is  greater  than  that  In  region  2.  Because  of 
this  difference  In  overpressure,  he  would  expect  that  the  Mach  shock 
would  rotate  forward. 

Laboratory  observations  have  shorn  that  the  Mach  stem  can  be 
considerably  curved  and  tends  to  bulge  forward.  This  Implies  that  the 
velocity  of  the  Mach  stem  Is  greater  than  In  the  Incident  shock  In  the 
direction  parallel  to  the  wedge  surface.  Thus  to  the  observer  on  the 
triple  point,  the  Mach  stem  should  toe  out  and  appear  to  rotate  or 
develop  curvature.  Both  effects  have  been  observed.  Sternberg  studied 
the  curvature  In  Mach  shocks  and  proposed  that  as  a  reason  for  the 
apparent  lack  of  agreement  between  theory  and  experiment  (61:182).  The 
towing  out  phenonmenon  can  be  seen  clearly  In  many  photographs  of  Mach 
stems  (39:flg  7-11, k,l)  (17:b43), 

The  actual  angle  of  Incidence  of  the  ambient  air  on  the  Mach  stem  Is 
more  correctly  described  by  an  angle  0-<T  where  O  represents  an  amount 
of  rotation  or  curvature.  However,  this  destroys  the  self-similarity 
assumption  made  earlier.  The  airflows  will  depend  on  the  degree  of 
rotation  about  the  triple  point  and  which  changes  with  time  and  distance 
traveled.  The  introduction  of  another  variable  requires  another 
condition  for  solution  of  the  system.  This  causes  the  same  problem  as  In 
the  solution  to  the  three  shock  system  before.  There  will  always  be  too 
many  variables  and  too  few  conditions  unless  some  assumption  Is  made.  It 
is  useful  to  recognize  at  this  point  that  the  Mach  stem  can  be  expected 


to  toe  out  and  not  remain  vertical  to  the  ground  in  the  region  of  the 
triple  point. 

From  the  frame  of  reference  of  an  observer  located  at  the  base  of  the 
Mach  stem  (point  E  in  figure  8),  the  flow  of  air  must  be  parallel  to  the 
wall  at  all  points.  If  the  shock  was  deflected  towards  the  ground,  a 
reflected  shock  would  be  generated  at  some  point  to  the  rear  of  the  Mach 
shock.  Similarly  If  the  flow  were  diverted  upwards,  It  would  cause  a  low 
pressure  area  to  be  generated  to  the  rear  of  the  Mach  shock.  The 
requirement  for  flow  parallel  to  the  wall  requires  that  the  Mach  stem  be 
perpendicular  to  the  wall  In  the  region  near  point  E.  Thus  regardless  of 
the  angle  of  the  Mach  stem  In  the  vicinity  of  the  triple  point,  the  Mach 
stem  must  Intersect  the  wall  surface  at  a  90  degree  angle.  Given  the 
fact  that  the  Mach  stem  rotates  about  the  triple  point,  the  Mach  stem 
must  have  a  curvature  that  changes  with  time. 

By  the  same  logic  as  that  which  defined  the  angle  in  equation  IV-5, 
or  may  be  defined  with  the  following  relation: 

cot  or  *  K3tan^  (IV-8) 

The  strength  of  the  reflected  shock  is  defined  as  ^  Because 

the  pressure  across  the  slipstream  is  constant,  ^3“?^ .  Thus  the  strength 
of  the  reflected  wave  can  be  expressed  in  terms  of  the  pressures  in 
regions  2  and  4  which  can  be  found  from  the  oblique  shock  relations. 


The  ratio  of  the  velocities  of  the  air  in  regions  1  and  2  are  given  by 
Mirels  as  (44:bl): 

I2  “  *  1  (IV-10) 

Ul/  Klsin/3 

Mirels  also  gives  a  relation  developed  using  the  method  of 
characteristics  for  the  flow  over  wedges  of  half-angle  A  in  free  flowing 
air.  This  equation  can  be  used  to  unite  the  preceedlng  relations. 

tan* A  =  rciC2M.,»  -  l)  ,  i  ^  +  (IV-11) 

)  1  +  CU  f  \  Mp-  _  y  +  1  )• 

Where  Cl  =  (  y  +  l)/(  y  -  1) 

C2  »  2  y/(  y  +  1) 

Using  the  above  relations,  equation  IV-11  can  be  expressed  as  a 
function  of  the  variables  B  and  ^  only.  Because  of  the  complexity  of 
the  relation,  it  has  not  been  possible  to  separate  the  equation  to  find 
^  directly.  Even  if  it  could  be  done,  the  parameters  K1  through  K4  are 
also  functions  of  the  angle  ^  .  Hence  some  iterative  method  of  solution 
must  be  used  to  find  this  angle.  The  relation  expressed  in  terms  of  0 
and  ^  is  as  follows: 

<Klcot0  +  KlK3tan^  -t-  cot(  g  +  VO [cot tan^  +K3]  )  * 

(  K1K3  +  Klcot0tan^  +  cot(9  +  ^HK3tan9  -  couf  ]  ) 

-  <C1C2(H-Kl*tan*(  ^  +  -  1>  f  KlK2fK2/K4  -  1) _  )* 

(  (1  +  C1K2/K4)K1K2  P)  [1  +  Kinan*(  ^  +  ^^)]-K2/K4+lf 

(IV-12) 


The  above  relation  is  used  to  find  the  angle  for  the  situation  of  a 


strong  shock,  ^  <  .15.  For  weaker  shocks  (  ^  >.15),  the  theory 
predicts  a  much  higher  angle  0  than  Is  experimentally  found  and  predicts 
that  regular  reflection  terminates  at  an  earlier  angle  of  incidence  than 
experments  show.  Thus  some  other  means  of  predicting  this  angle  is 
required  in  the  region  of  weak  shocks. 

Use  of  Experimental  Data. 

Experimental  data  was  used  by  this  author  to  derive  a  set  of  curves 
which  would  be  used  to  predict  the  angle  ^  given  an  incident  shock 
strength  and  angle  of  Incidence.  The  data  used  was  a  set  of  fitted 
curves  for  0  given  by  Harlow  (29:56)  and  supplemented  by  a  series  of 
experiments  done  at  the  University  of  Toronto  (5:table  5),  (17:bl)  in  air 
and  nitrogen.  Curves  were  fitted  to  the  following  shock  strengths,  .15, 
.30,  .40,  .50,  .60,  .70,  .80,  '=nd  .90.  Shocks  with  strengths  between  the 
fitted  values  are  found  by  linear  interpolation. 

The  set  of  curves  was  checked  against  the  actual  data  for  accuracy  and 
found  to  predict  the  value  of  to  within  ten  percent.  The  comparison 
in  accuracy  is  particularly  striking  with  respect  to  the  prediction  of 
the  onset  of  Mach  reflection.  For  example,  the  computational  technique 
presented  above  predicts  Mach  reflection  to  begin  at  a  Mach  number  of 
about  1.325  for  a  wedge  angle  of  40  degrees.  This  model  predicts  the 
onset  at  about  1.175.  The  experimental  value  is  a  Mach  number  of 
approximately  1.2.  At  higher  Mach  numbers  the  two  predictions  merge  to 
the  same  answer. 

The  set  of  curves  used  here  could  be  substantially  Improved  if 


additional  detailed  data  was  available  for  weak  incident  shocks.  The 
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Toronto  experiments  covered  wedge  angles  of  two  degrees  to  50  degrees  and 
Mach  numbers  ranging  from  1.6  to  6.9,  Referring  to  figure  19,  it  can  be 
seen  that  at  low  Mach  numbers,  small  errors  can  make  much  larger  errors 
in  the  value  of  Because  of  the  way  the  data  is  used,  more  detailed 

Information  from  experiments  at  low  Mach  numbers  and  small  wedge  angles 
would  enhance  the  predictions  considerably. 


Solution  for  Heights  Greater  than  Height  of  Burst, 


when  the  height  of  the  triple  point  reaches  the  height  of  burst,  the 
shocks  assume  a  new  configuration  (see  figure  12).  This  configuration  is 
not  reproducible  in  the  laboratory  with  planar  shocks  on  wedges.  The 
Mach  shock  theory  described  before  does  not  work  well  for  angles  of 
incidence  greater  than  85  degrees  (40:3,8).  The  equations  do  not  give 
meaningful  answers  for  shocks  with  angles  of  Incidence  greater  than  90 
degrees  either.  Thus  a  new  method  of  examining  the  problem  was  developed 
by  this  author. 

Generally,  certain  assumptions  are  made  to  make  the  problem  solvable. 
Both  the  incident  and  the  reflected  shocks  are  assumed  to  be  spherical  in 
shape.  The  inhomogeneity  of  the  atmosphere  and  its  effect  on  the  shape 
of  the  incident  shock  are  Ignored.  The  large  variations  in  pressure  and 
density  that  occur  in  the  region  behind  the  incident  shock  may  cause 
large  deviations  from  spherical  sjrmmetry  in  the  reflected  shock, 
especially  near  ground  zero.  However,  the  triple  point  reaches  the 
height  of  burst  at  a  ground  distance  greater  than  five  to  six  times  the 
height  of  burst.  By  this  time,  the  Incident  shock  has  lost  much  of  its 
strength  and  Is  barely  sonic  for  scaled  heights  of  burst  greater  than  200 


Figure  12.  Shock  Configuration  at  Late  Times 
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feet.  Because  of  this  fact,  the  shock  can  be  approximated  as  spherical 
and  an  acoustic  reflection  approximation  made.  For  weak  shocks,  the 
angle  of  reflection  is  equal  to  the  angle  of  incidence.  Similarly,  the 
overpressure  in  the  reflected  shock  is  equal  to  the  incident  shock.  From 
the  acoustic  approximation,  the  shape  and  strength  of  the  reflected  shock 
can  be  found  from  an  'image  burst'  fired  at  a  point  one  height  of  burst 
below  the  surface  of  the  earth  at  ground  zero.  The  overpressure  of  the 
reflected  shock  is  assumed  to  be  the  same  as  a  shock  traveling  in  free 
air  from  the  image  burst.  Simple  geometrical  relations  may  then  be  used 
to  find  the  angle  of  Incidence  and  the  angle  between  the  two  shocks. 

Once  these  angles  are  known,  the  shocks  are  treated  as  planar  shocks  for 
that  instant  in  time,  and  the  angle  calculated. 

The  movement  of  the  Intersection  of  the  two  planar  shocks  is 
illustrated  in  figure  13.  The  velocities  of  the  incident  and  reflected 
shocks  are  designated  Vj^  and  V2  respectively.  The  x  and  y  directions  are 
parallel  and  perpendicular  to  the  plane  of  the  earth's  surface.  The 
amount  of  movement  in  each  of  these  directions  during  an  infinitesimal 
time  dt  can  be  found  from  the  following  relations: 

(IV-13) 
(lV-14) 

The  velocity  of  each  shock  can  be  measured  easily  in  terms  of  the  Mach 
number  of  each  shock  which  in  turn  can  be  found  from  the  overpressure. 
However,  the  reflected  shock  travels  faster  in  the  air  behind  the 


dx  =  COS0  V^dt  -  cos(0+  6  )  V2dt 
sin  i  sin  (J 

dy  =  sin(0  +  6  )  V^dt  -  sin0  V^dt 
sin  i  sin  0 


Incident  shock  than  it  could  in  ambient  air.  Thus  the  true  velocity  of 
the  reflected  shock  is  not  given  by  the  image  burst.  The  difference  in 
velocity  can  be  estimated  from  the  ratio  of  the  sound  speeds  in  the 
shocked  and  unshocked  air.  Using  the  Rankine-Hugoniot  relations,  the 
following  relationship  was  found: 

/l^f  =  (Po+dP)(7Po+dP)  =  f^  (IV-15) 

\Sjy  Po(7Po+6dP) 

Where:  Po  =  ambient  atmospheric  pressure 

dP  =  overpressure  in  the  incident  shock 
S2  =  sound  speed  in  the  shocked  air 

=  sound  speed  in  the  ambient  atmosphere 

It  should  be  noted  .that  the  speed  of  sound  in  the  shocked  air  may  vary 
greatly  from  area  to  area.  The  ratio  above  is  an  attempt  to  quantify  the 
velocity  of  the  reflected  shock  in  the  area  immediately  behind  the 
incident  shock.  Using  the  above  relation  and  the  fact  that  tan ^  = 
dy/dx,  a  relationship  may  be  written: 

atn  rM,fsin(^+j)  -  Misin^l  (IV-16) 

L^icos^  -  M2fcos(0+  i  )J 

Where:  (0+i)  =  180  -  angle  of  Incidence 

0 =  90  -  atn( (HTP+HOB) /ground  distance] 
f  =  ratio  of  sound  speeds 

The  relations  immediately  preceeding  may  be  used  whenever  the  height 
of  the  triple  point  is  greater  than  the  height  of  burst.  The  angles 
generated  by  this  method  compare  favorably  to  the  angles  found  by 
anylysls  of  the  triple  point  path  data  in  the  Reflect-4  computer  codes 


(58:147-173).  Reflect-4  Is  a  two-dimensional  hydrodynamics  code  for  the 
analysis  of  Mach  reflection.  The  angles  found  by  the  method  above  are 
consistently  less  than  the  reflect  code  would  Indicate  by  about  20 
percent.  This  would  seem  to  Indicate  that  the  reflected  shock  wave 
travels  faster  than  the  estimate  based  upon  the  ratio  of  sound  speeds. 

In  evaluating  the  above  set  of  equations,  there  are  two  limiting  cases 
that  they  must  converge  to.  For  low  heights  of  burst,  the  triple  points 
on  each  side  of  the  burst  will  rise  very  quickly  and  meet  above  the 
burst,  creating  a  shock  system  similar  to  a  contact  surface  burst.  For 
very  high  heights  of  burst  and  at  long  distances,  the  reflected  shock 
will  be  very  weak  in  relation  to  the  Incident  shock.  Thus  the  angle  ^ 
will  converge  to  some  value  and  not  change  significantly.  The  method 
described  here  does  predict  that  for  low  heights  of  burst,  the  triple 
points  will  arch  back  and  meet  above  the  burst.  However,  In  the  case  of 
high  heights  of  burst  and  long  distances,  the  equations  predict  that  the 
angle  0  will  continue  to  grow  in  magnitude.  Thus  the  method  described 
above  should  be  used  with  care  since  It  appears  not  to  be  universally 
valid  for  all  situations. 


Comparison  to  Nuclear  Explosions 

The  computer  program  In  Appendix  F  utilizes  the  methods  In  this 
chapter  to  trace  the  trajectory  of  the  triple  point  in  space.  The 
equations  given  are  used  In  each  of  the  three  regions  to  predict  the 
angle  0  .  This  angle  Is  then  used  In  equation  IV-1  and  Integrated  using 
Simpson's  rule.  Two  scaled  heights  of  burst  were  selected  for 
comparison.  Teapot  number  4  and  number  12  were  two  shots  with  a  135  foot 
scaled  HOB.  The  path  of  the  triple  point  was  determined  In  detail  for 


these  two  shots  and  present  a  good  example  of  the  degree  of  correlation 
between  any  estimates  and  an  actual  detonation.  The  second  comparison  is 
made  for  a  700  foot  scaled  HOB  using  the  Reflect-4  Code  determinations. 
Reflect-4  is  a  two-dimensional  hydrodynamic  code  used  for  Mach  stem 
studies  by  the  Kaman  Avidyne  Corporation.  The  comparison  with  this 
prediction  will  enable  a  comparison  with  current  state  of  the  art  in 
shock  modeling.  In  each  case,  the  prediction  based  upon  Potocki's 
curve-fits  to  EM-1  will  be  used  to  show  how  well  current  practice  relates 
to  the  same  examples.  The  curves  in  EM-1  are  the  basis  for  most  of 
today's  predictions  of  the  triple  point  trajectory.  These  curves  are 
directly  traceable  to  the  work  of  Hesse  and  Kelso  in  1955  (35:11,22,23). 

Figure  14  shows  the  comparison  between  the  two  Teapot  detonations  and 
the  program  predictions  for  a  135  foot  scaled  height  of  burst.  The 
prediction  made  by  the  program  is  too  low  at  early  times  and  too  high  at 
much  later  times.  To  some  extent,  the  rapid  rise  of  the  triple  point 
from  the  Teapot  detonations  is  due  to  precursor  activity.  The  Mach  stem 
formed  much  earlier  than  predicted  and  rose  much  faster.  This  may 
account  for  some  of  the  discrepancy.  There  is  a  large  percent  difference 
in  triple  point  height  just  between  the  asphalt  and  soil  surfaces  in  the 
Teapot  12  tests.  The  water  surface  on  that  test  would  correspond  to  the 
most  ideal  surface  expected.  The  calculations  are  much  lower  than  even 
this  surface. 

At  a  700  foot  scaled  height  of  burst  the  prediction  does  much  better 
(Figure  15).  In  general,  the  prediction  is  better  than  Potocki's 
prediction.  The  change  in  slope  of  the  prediction  at  the  altitude  of  the 
burst  is  caused  by  the  change  in  routines  from  experimental  data  to  the 
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expression  used  to  find  ^  at  altitudes  above  the  triple  point.  The 
experimental  data  used  is  weak  in  this  area  and  needs  to  be  improved. 


The  results  for  the  135  foot  scaled  height  of  burst  are  indicative  of 
the  problem  of  prediction  of  the  triple  point  trajectory.  For  the  same 
burst,  but  for  different  ground  types,  the  path  was  very  different. 

There  is  also  the  large  difference  between  Teapot  4  and  Teapot  12.  They 
were  fired  in  similar  atmospheric  conditions,  over  the  same  type  of  soil 
and  in  the  same  season  of  the  year.  The  actual  yield  of  the  weapon  and 
height  of  burst  were  different  and  caused  a  large  difference  in  the 
triple  point  heights  at  the  same  scaled  range.  It  seems  reasonable  to 
conclude  that  errors  on  the  order  of  ten  to  twenty  percent  in  the  height 
of  the  triple  point  are  reasonable  at  the  level  of  prediction  ability 
available  today.  The  addition  of  more  physics  into  the  prediction,  and 
taking  into  account  the  precursor  heating  of  the  air,  the  change  In  air 
density  with  temperature  and  altitude  may  enable  a  better  prediction  of 
the  path.  However,  the  addition  of  these  effects  means  that  the 
prediction  will  take  considerably  more  time  to  calculate  and  will  be 
device  and  location  dependent.  These  factors  are  not  desirable  for  many 
uses . 
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L  Empirical  Results  for  Shape.  Ta,  and  Overpressure. 


To  the  military  targeteer,  the  effects  that  must  be  known  in  order  to 
predict  the  damage  to  aircraft  and  missiles  in  the  path  of  the  Mach  stem 
include  the  shape  or  orientation  of  the  shock  wave,  the  overpressure,  the 
dynamic  pressure,  and  the  time  of  arrival.  These  effects  when  combined 
with  the  type  and  orientation  of  the  aircraft  allow  determination  of  the 
probability  of  survival.  Presently  the  Mach  stem  is  assumed  to  be 
vertical  at  all  points  above  the  ground.  The  time  of  arrival  and 
overpressure  may  be  approximated  by  the  time  of  arrival  and  overpressure 
on  the  ground.  The  overpressure  at  altitude  is  scaled  using  Sachs 
scaling.  Current  shock  theory  is  not  capable  of  predicting  any  of  these 
effects  at  any  point  for  any  geometry. 

Shape  and  Orientation  of  the  Shock  Wave. 

The  examination  of  photographs  of  above-ground  nuclear  bursts  as  well 
as  data  from  hydrodynamic  calculations  such  as  Reflect-4  clearly 
demonstrate  that  the  assumption  of  a  vertical  Mach  shock  is  wrong. 
Examination  of  the  data  shows  that  there  are  in  fact  three  distinct 
shapes  to  the  Mach  stem  (See  figure  16).  At  early  times,  the  triple 
point  leads  the  base  of  the  Mach  stem  by  as  much  as  ten  percent  of  the 
ground  range.  For  low  heights  of  burst  such  as  200  scaled  feet,  the 
triple  point  may  lead  as  much  as  fifteen  to  seventeen  percent 
(59:30,150).  During  this  period,  the  Mach  shock  is  concave  Inward.  As 
the  triple  point  approaches  the  height  of  burst,  the  base  of  the  Mach 
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shock  catches  up  and  the  Mach  stem  Is  approximately  vertical.  As  the 
height  of  the  triple  point  continues  to  rise,  the  base  of  the  Mach  stem 
begins  to  lead  the  triple  point.  As  time  goes  on,  the  shape  of  the 
entire  shock  system  of  the  Incident  shock  and  Mach  shock,  approaches  a 
hemispherical  shape  centered  at  ground  zero.  Because  of  the  Increase  In 
speed  of  the  shock  at  low  heights  of  bursts  due  to  the  layer  of  heated 
air  near  the  ground,  this  hemispherical  shape  can  be  Introduced  by  the 
time  the  Mach  shock  Is  50  feet  high  (45:166,201).  In  other  words,  the 
presence  of  precursor  conditions  may  enhance  the  tendency  towards  the 
hemispheric  shape. 

The  airflow  In  the  vicinity  of  the  triple  point  requires  that  the  mach 
stem  lead  the  Incident  shock  when  the  triple  point  Is  above  the  height  of 
burst.  In  this  case,  the  Incoming  air  Is  deflected  upward  by  the 
Incident  shock  (See  figure  17).  The  upward  deflection  Is  caused  by  the 
angle  of  Intersection  between  the  Incoming  air  and  the  shock  (Designated 
0  In  figure  17).  Where  before  the  angle  was  acute,  the  angle  In  the  new 
configuration  Is  obtuse.  The  reflected  shock  also  has  an  obtuse 
Intersection  angle  and  deflects  the  flow  of  air  upward.  If  the  mach 
shock  Is  to  deflect  the  air  through  the  same  angle  as  the  other  two 
shocks  (1.  e.  ^2)*  It  must  lead  the  Incident  shock.  If  the  shock 

was  vertical,  the  air  would  be  deflected  down  by  the  mach  shock,  causing 
the  airflow  behind  the  triple  point  to  be  discontinuous.  This 
discontinuity  does  not  occur.  As  the  three  shocks  lose  their  energy,  the 
reflected  shock  becomes  weakest  first.  When  It  is  very  weak,  the 
deflection  angle  of  the  air  due  to  It  will  become  negligible.  When  that 
happens,  the  deflection  caused  by  the  incident  shock  must  equal  the 


deflection  of  the  mach  shock.  At  this  time  the  tangents  of  the  two 
shocks  must  converge. 


Figure  17.  Airflows  Near  Triple  Point  When  HTP  is  Greater  than  HOB 


Since  damage  due  to  overpressure  for  aircraft  occurs  in  the  1  to  20 
psi.  overpressure  range,  the  region  of  most  interest  is  this  distant 
region.  The  equations  given  here  are  valid  when  the  triple  point  is 
heigher  than  the  height  of  burst.  The  effective  radius  of  the  Mach  stem 
at  these  overpressures  is  given  in  table  1. 


Scaled 

Scaled 

Scaled  HOB 

Ground  Range 

Effective  Radius 

Ratio 

(feet) 

(feet) 

(feet) 

200 

705 

728 

1.1 

400 

4,982 

9,878 

1.98 

700 

5,614 

15,833 

2.77 

Data  from  analysis  of  reference  59. 


Table  1.  Effective  Radius  of  Mach  Stem 
The  radius  of  the  Mach  stem  decreases  more  slowly  at  higher  scaled 


heights  of  burst  because  it  grows  in  height  many  times  more  slowly.  The 
lower  the  height  of  burst,  the  sooner  it  resembles  a  contact  surface 
burst  to  a  distant  observer.  In  fact,  the  Reflect  codes  predict  that  for 
a  200  foot  scaled  height  of  burst,  the  triple  point  will  continue  to  rise 
after  surpassing  the  height  of  burst  elevation,  but  will  decrease  in 
ground  range.  Then  the  entire  shock  front  rapidly  becomes  reinforced  so 
that  within  one-half  second  of  the  triple  point  passing  the  height  of 
burst,  the  Mach  stems  meet  at  the  apex  of  the  shock  front  (59:151).  The 
effective  radius  of  the  Mach  stem  at  large  distances  from  ground  zero  may 
be  approximated  by  the  following  empirical  formula,  which  is  a  curve  fit 
to  the  data  in  Table  1. 

Eff.  Radius  =  Re  =  ( .0000573SHOB^+.317SHOB+.574)(Gnd.  Range)  (V-1) 

Where  SHOB  is  Scaled  HOB,  measured  in  hundreds  of  feet. 

Once  the  radius  is  known,  the  orientation  of  the  Mach  stem  can  be 
found.  From  this,  the  direction  of  the  impulse  imparted  by  the  shock 
front  with  respect  to  the  ground  may  be  calculated.  This  direction  is 
also  the  direction  of  the  'gust*  that  is  used  in  the  calculation  of  the 
wing  loading.  T  as  given  in  the  following  equation  is  the  direction 
normal  to  the  shock  surface. 

T  =  arctan/Helght  of  Target  \  (V-2) 

\Effective  Radius/ 

Time  of  Arrival. 

In  attempting  to  find  a  relation  that  predicts  the  time  of  arrival, 


two  basic  relations  were  used.  Both  Horizons  Technology  (51:82)  and 
Brode  (8:15)  give  general  relations  for  time  of  arrival  in  free  air.  The 


Horizons  Technology  equation  is  essentially  a  curve  fit  to  the  one 


kiloton  standard.  To  properly  scale  back  to  the  desired  case,  detailed 
knowledge  of  the  temperature  and  pressure  of  the  target  is  required. 

Even  with  this  knowledge  for  the  original  bursts,  this  equation  failed  to 
predict  the  actual  time  of  arrival  as  well  as  the  equation  given  by  Erode 
(See  table  2).  The  best  fit  found  to  predict  the  actual  times  of  arrival 
is  given  by  Erode:  (8:15) 


Ta  =  (. 54291-21. 185R+361.81R^+2383R^)W^^^  msec  (V-3) 

(1  +  2.04797R  +  2.68717R^) 

Where:  R  =  Slant  range  to  target  (kilofeet)  scaled  by  the  yield 

and  pressure  on  the  ground. 

W  =  Yield  (kilotons) 

Ta  =  Time  of  arrival  (milliseconds) 

This  formula  is  accurate  to  within  7.7  percent  of  the  actual  values. 
Since  the  original  slant  ranges  are  known  to  within  eleven  percent, 
greater  accuracy  cannot  be  expected.  In  all  cases  but  one,  the  formula 
predicts  a  slightly  earlier  time  of  arrival  than  is  observed.  The  sole 
exception  is  the  Ivy  King  test  in  which  the  actual  time  of  arrival  is 
consistently  earlier  than  predicted.  This  burst  differed  from  the  others 
in  that  it  was  550  Kt .  and  the  others  were  43  Kt.  or  less.  The  Ivy  King 
test  also  differs  from  the  other  tests  in  that  it  appeared  to  have  a  much 
greater  blast  efficiency  than  previous  tests.  The  free  air  and  Mach  stem 
overpressures  were  much  greater  than  expected.  The  greater  overpressure 


would  cause  a  decrease  in  time  of  arrival.  Because  cannister  data  is 
available  for  only  one  high  yield  weapon,  it  is  not  possible  to  tell  if 
the  Ivy  King  data  is  a  device  peculiar  effect  or  to  be  expected  of  all 
high  yield  weapons. 


Author 

Yield 

Min  Error 

Max  Error 

Ave  Error 

Erode 

W 

.22  % 

7.69% 

4.67  % 

Erode 

2W 

.377% 

14.3  % 

8.53  % 

Horiz  Tech 

W 

.23  % 

34.5  % 

13.58  % 

Horiz  Tech 

2W 

1.3  % 

33.5  % 

11.16  % 

Table  2.  Comparison  of  Time  of  Arrival  Equations  for  Mach  Stem  Data 

The  Mach  shock  seems  to  propagate  at  approximately  the  same  speed  as  a 
free  air  shock  at  long  distances.  As  shown  earlier,  at  long  distances 
the  Mach  shock  assumes  the  hemispherical  shape  of  the  incident  shock. 

Thus  the  entire  shock  system  appears  to  a  far  observer  to  have  a 
hemispherical  shape.  Under  such  conditions,  the  time  of  arrival  of  the 
Mach  shock  would  seem  to  be  the  same  as  a  free  air  shock  at  the  same 
slant  range.  The  equation  given  is  valid  then  for  arrival  times  greater 
than  three  seconds  or  when  the  height  of  the  triple  point  is  much  greater 
than  the  height  of  burst. 


Overpressures . 

There  are  only  two  points  in  the  Mach  stem  at  which  the  overpressure 
may  be  estimated  with  any  degree  of  certainty.  These  are  the 


overpressure  at  the  base  of  the  Mach  stem  and  the  overpressure  just  below 
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the  triple  point.  The  overpressure  at  the  base  of  the  Mach  stem  is  known 
from  ground  measurements  made  during  the  above  ground  nuclear  tests.  The 
overpressure  may  also  be  estimated  by  semi-empirical  methods  based  upon 
the  free  air  overpressure  at  equal  ranges  with  considerable  accuracy 
(8:14-27).  However,  in  almost  all  cases,  the  overpressure  measurements 
were  made  at  real  ground  ranges  of  less  than  5000  feet.  The  height  of 
burst  curves  in  EM-1  are  accepted  as  the  standard  for  the  estimation  of 
the  ground  overpressure.  They  also  are  limited  to  a  scaled  ground  range 
of  4000  feet.  Mathematical  curve  fits  made  from  this  set  of  curves  may 
be  used  to  extend  the  predictions  with  an  accuracy  of  fifteen  percent 
(See  for  example  (52:42)).  At  the  triple  point,  the  overpressure  may  be 
found  from  equation  IV-7.  Since  the  Mach  stem  is  not  vertical,  the  true 
angle  of  incidence  is  not  90- (p  ,  but  some  angle  90-lp-ffy  where  <T  is 
the  Mach  stem's  deviation  from  vertical.  Thus  for  any  given  instant  in 
time,  one  may  find  the  value  of  the  overpressure  just  below  the  triple 
point  if  he  knows  rate  of  rise  of  the  triple  point  and  the  amount  of 
deviation  from  vertical  in  the  Mach  stem. 

However,  the  military  planner  needs  to  know  the  variation  in  pressure 
with  height.  The  fact  that  the  Mach  stem  changes  curvature  in  itself  is 
an  indication  that  there  is  a  pressure  distribution  along  the  Mach  front. 
Analytical  theory  that  would  predict  the  physical  characteristics  of  the 
Mach  shock  is  nearly  nonexistent  in  planar  shock  theory  and  completely 
nonexistent  in  the  case  of  spherical  shocks.  In  1951  Fletcher,  Bleakney, 
and  Taub  reviewed  three  theories  that  purported  to  describe  conditions 
within  planar  shocks  (26:271).  Each  of  these  theories  was  valid  only  in 
the  region  of  vanishingly  weak  shocks  and  angles  of  incidence  of  nearly 


90  degrees.  Since  that  time,  there  have  been  no  significant 
improvements.  The  most  promising  of  the  theories  was  one  by  Bargmann  who 
predicted  that  the  pressure  in  the  Mach  stem  decreases  monotonically  with 
increasing  vertical  height  along  the  shock  and  specified  the  shape  of  the 
pressure  contours  behind  the  Mach  stem  (4:43).  The  characteristic  shape 
of  these  contours  has  been  experimentally  verified  and  taken  as  evidence 
for  his  theory.  Apart  from  the  expectation  of  a  pressure  variation  with 
altitude,  this  theory  is  of  little  use  over  the  wide  range  of  shock 
strengths  and  angles  of  incidence  to  be  found  in  a  nuclear  detonation. 

During  the  atmospheric  nuclear  testing  in  the  1950s,  cannisters 
occasionally  were  dropped  to  measure  overpressures  at  altitude.  Later, 
more  sophisticated  photographic  techniques  were  used  to  measure  free  air 
shock  velocities  and  overpressures.  A  total  of  40  unclassified  cannlster 
measurements  were  made  in  the  Mach  region.  Two  more  measurements  are 
still  classified.  When  the  measurements  were  compared  to  the  standard 
method  of  prediction,  the  measurements  differed  from  the  predictions  in  a 
consistent  manner.  Figure  18  is  a  plot  of  the  actual  data  in  terms  of 
the  fractional  distance  of  the  gage  height  up  the  Mach  stem  versus  the 
measured  overpressure  divided  by  the  predicted  overpressure.  The 
predicted  overpressure  was  found  by  using  the  equations  given  by  Horizons 
Technology  (51:74-5).  The  height  of  burst  in  each  case  was  scaled  using 
the  designated  yield  and  ambient  air  pressure  on  the  surface.  The 
equations  gave  the  overpressure  on  the  ground  for  the  1  klloton  reference 
case.  This  pressure  was  then  scaled  to  the  altitude  of  the  particular 
measurement  using  Sach's  scaling.  The  ambient  air  pressures  on  the 
ground  and  at  altitude  are  known  in  all  cases  from  radiosonde  readings 


fig.  18  Variation  of  Pressure  up  the  Mach  Stem 
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before  the  tests.  If  the  data  agreed  with  the  predictions,  it  would 
follow  the  line  y=l.  However,  at  heights  of  about  20  percent  of  the 
triple  point  height,  the  overpressure  is  significantly  greater  than 
expected.  At  locations  approaching  the  triple  point  height,  the 
overpressure  drops  off  much  faster  than  expected.  As  the  height  of  the 
triple  point  and  ground  range  increases,  the  overpressure  near  the  triple 
point  can  be  expected  to  drop  off  because  of  the  weakening  of  the 
incident  shock  and  the  very  great  angle  of  incidence. 

Appendix  B  is  a  summary  of  data  that  supports  the  finding  of  the 
described  pressure  distribution.  This  effect  of  a  pressure  increase 
above  the  ground  followed  by  a  decrease  is  found  not  only  in  nuclear 
data,  but  also  in  small  lb.  TNT)  explosions  and  larger  (1000 
lb.  Penolite)  explosions.  The  effect  is  also  seen  in  some  of  the  planar 
shock  data  that  was  done  at  the  University  of  Toronto.  The  fact  that 
this  effect  is  found  in  Mach  shocks  other  than  nuclear  allows  several 
possible  causes  to  be  eliminated.  The  small  HE  tests  Indicate  that  the 
effect  is  not  caused  by  atmospheric  inversions,  funneling,  or  defects  in 
the  Sachs  scaling  laws.  Between  the  four  types  of  experiments,  a  large 
variation  in  scale  is  covered.  The  shock  tube  Mach  stems  may  be  measured 
in  fractions  of  an  inch,  while  the  HE  tests  cover  a  span  of  several  feet 
and  the  nuclear  data  covers  a  distance  of  up  to  several  miles.  The  wide 
range  of  scales  would  seem  to  indicate  that  the  cause  of  this  effect  is 
not  due  to  atmospheric  irregularities  of  any  sort.  The  cause  must  lie  in 
the  way  in  which  the  Mach  stem  is  formed  and  the  resulting  pressure 
changes  that  occur  throughout  the  life  of  the  phenomenon. 

The  pressure  behind  the  Mach  stem  may  be  found  from  the  following 
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relation; 


Pj  =  Pj  2yM,  ^cos^0  -  y-  1  )  (V-4) 

((>'+  I)cos''(  ^  y+  1) 

Where  Instantaneous  angle  of  rise  of  the  triple  point 

B  =  Wedge  angle 

B  =  90  -  angle  of  incidence  for  spherical  shocks 
Pj^  =  Ambient  air  pressure 
^2  ~  Pressure  behind  the  Mach  stem 

An  increase  in  the  angle  ^  will  always  correspond  to  an  Increase  in 
the  pressure  behind  the  Mach  shock.  Similarly,  when  decreases,  the 
pressure  behind  the  Mach  stem  will  decrease.  An  increase  in  B  will  also 
cause  an  Increase  in  the  pressure  if  all  other  factors  are  held  constant. 

Figure  19  is  a  graph  showing  how  the  angle  ij/  varies  experimentally 
with  the  wedge  angle  and  Mach  number.  Because  the  angle  tlf  governs  the 
angle  at  which  the  incoming  air  strikes  the  Mach  stem,  an  increase  in 
translates  into  an  Increase  in  overpressure  in  the  Mach  stem  at  that 
point.  Thus  the  vertical  axis  is  also  related  to  overpressure  in  the 
Mach  stem. 


TIME 


Figure  20.  Pressure  Increase  in  the  Mach  Stem  vs.  Time 
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In  a  spherical  shock,  the  angle  B  will  decrease  with  time  due  to  the 
rise  of  the  Mach  stem.  As  shown  in  figure  20,  this  will  cause  a  gradual 
decrease  in  the  pressure  in  the  mach  stem  just  below  the  triple  point. 

The  angle  0  will  increase  with  time  because  the  Mach  number  of  the  shock 
decreases  with  time  and  because  the  angle  B  is  decreasing.  This  will 
cause  the  overpressure  just  below  the  triple  point  to  rise  (See  figure 
19).  The  net  result  of  the  two  opposing  forces  is  the  third  curve  in 
figure  20.  The  overpressure  in  the  mach  stem  in  relation  to  the 
overpressure  in  the  incident  shock  first  Increases  and  then  decreases. 
Thus  one  would  expect  a  'bulge*  in  the  pressure  distribution  in  the  Mach 
stem.  The  net  effect  is  to  have  some  overpressure  on  the  ground,  an 
Increase  in  pressure  at  some  height  above  the  surface,  and  above  that  a 
decrease  in  pressure. 

At  this  point  it  must  be  noted  that  the  time  of  measurement  is 
critical  to  the  measurement  of  the  pressure  distribution.  At  very  early 
times,  only  a  pressure  increase  would  be  measured.  At  very  late  times, 
the  length  of  the  decrease  in  pressure  will  be  large  in  proportion  to  the 
length  of  the  Mach  stem.  At  middle  times,  the  region  of  Increase  in 
pressure  may  be  proportionally  large  and  gradually  diminish  in  size.  It 
is  also  clear  that  the  point  at  which  the  maximum  overpressure  occurs  in 
terms  of  the  fractional  distance  up  the  Mach  stem  changes  in  time.  At 
the  ground  distance  at  which  a  military  planner  would  require  the 
overpressure,  the  location  of  this  maximum  is  not  likely  to  change 
radically  with  time. 

A  nuclear  burst  is  characterized  by  a  spherical  symmetry  which  causes 
the  shock  to  lose  energy  relatively  rapidly  and  also  allows  the  angle  of 
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Incidence  to  change  with  time.  Examining  the  graph  of  figure  19,  It  can 
be  seen  that  at  large  wedge  angles,  the  amount  of  Increase  In  0  Is 
small.  Similarly,  the  amount  of  Increase  In  small  wedge  angles  Is  very 
great.  Thus  In  a  low  height  of  burst  (200  scaled  feet  or  less)  where  the 
triple  point  rises  very  rapidly  and  thus  the  angle  of  Incidence  rises 
very  rapidly  (angle  of  Incidence  equals  90  -  wedge  angle).  It  Is  possible 
to  get  a  much  higher  Increase  In  pressure  than  In  a  higher  height  of 
burst. 

With  the  aid  of  figure  20,  It  can  be  shown  that  when  the  Mach  stem  Is 

generated  these  variations  of  pressure  exist.  With  time  these  pressure 

variations  will  tend  to  diffuse  through  the  shock  and  equalize  the 

overpressure  at  all  points  along  the  shock  front.  Whltham  has  shown  that 

the  difference  In  velocity  between  two  areas  of  different  overpressure 
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will  tend  to  decrease  as  t  (67:42-47,307-309).  Because  of  the  way 
overpressure  Is  related  to  shock  velocity  (See  Appendix  A),  the  resulting 
rate  of  decrease  of  the  difference  in  overpressure  Is  somewhat  slower. 

If  the  magnitude  of  the  overpressure  bulge  could  be  predicted,  the  time 
It  takes  to  diffuse  away  could  be  predicted  using  the  Whltham  theory. 
However,  at  this  time  the  magnitude  of  the  overpressure  bulge  is  not  able 
to  be  predicted. 

From  the  nuclear  data,  a  curve  fit  was  made  by  the  author  to  the  data 
In  figure  18.  This  represents  the  overpressure  profile  In  the  shock 
after  Sachs  scaling  has  been  applied.  This  curve  Is  useful  In  the  region 
where  the  profile  does  not  change  significantly.  It  applies  only  when 
the  height  of  the  triple  point  has  risen  to  at  least  twice  the  height  of 
burst.  This  region  corresponds  to  overpressures  of  20  psl.  or  less.  The 


curve  was  fit  by  a  sixth  order  polynomial  which  resulted  in  the  following 


equation: 


F  =  .993766  +  4.72762x  -  15.5804x*  -  21.66193x^  (V-6) 

+127.5835x^  -  155.5353X®  +  59.9366x® 


Where  x  =  Fractional  distance  up  the  Mach  stem;  Target  height/HTP 
F  =  Overpressure  ratio 


The  above  relation  may  be  used  to  calculate  the  overpressure  at  any 
given  point  in  the  Mach  stem  by  the  following  relations: 


Overpressure  at  altitude  =  (ground  overpressure)F(Pa)/Po 


Where  F  =  Overpressure  ratio  defined  above 


(V-7) 


Pa  =  Ambient  pressure  at  altitude 
Po  =  Ambient  ground  pressure 


The  overpressure  on  the  ground  may  be  found  from  any  of  several 
relations  including  the  graphs  in  EM-1,  Horizons  Technology  (51:73-75) 
and  by  Erode  (8:14-27). 

It  must  be  noted  that  the  overpressure  distribution  as  given  by 
equation  V-6  is  not  necessarily  the  overpressure  distribution  in  the  Mach 
stem  at  any  given  Instant.  It  actually  represents  the  variation  in 
pressure  that  would  be  seen  by  a  vertical  array  of  gages  as  the  shock 
wave  passed  by.  Because  of  the  hemispherical  shape  of  the  shock  wave, 
the  gages  on  the  bottom  of  the  array  would  be  struck  first  and  the 
topmost  ones  last.  Because  the  shock  wave  decays  in  the  small  amount  of 


time  between  the  striking  of  the  lowest  and  highest  gages,  the  actual 
overpressure  difference  In  the  Mach  stem  between  the  top  and  the  bottom 
may  not  be  as  large  as  the  distribution  Indicates. 

In  the  case  of  nuclear  detonations,  there  Is  an  additional  phenomenon 
which  occurs  that  may  reinforce  the  rise  In  overpressure  above  the 
surface.  The  radiation  from  a  nuclear  device  creates  a  layer  of  heated 
air  on  the  surface  of  the  earth.  This  situation  Is  similar  to  and  Is 
modeled  by  a  two  gas  system  (See  for  example  ref. 36).  A  light  gas  such 
as  helium  or  freon  Is  deposited  In  such  a  way  that  two  layers  are  formed, 
one  layer  less  dense  than  the  other.  When  a  shock  passes  through  this 
system,  there  will  be  a  large  Increase  in  overpressure  at  the  boundary 
between  the  two  layers.  It  can  be  shown  that  this  pressure  Increase 
occurs  In  nuclear  detonations  and  that  the  Increase  In  overpressure  above 
the  ground  extends  far  beyond  the  region  of  the  precursor  (See  Appendix  D 
for  detailed  data  from  several  nuclear  burets). 

The  results  of  this  effect  on  the  overpressure  distribution  described 

above  Is  twofold.  First,  the  tendency  for  the  overpressure  to  rise  above 

the  ground  Is  reinforced.  In  the  case  of  a  relatively  high  height  of 

burst,  the  overpressure  due  to  this  effect  will  be  greater  than  that  due 

to  the  change  in  the  angle  ^  .  Secondly,  this  Increase  In  overpressure 

tends  to  minimize  the  diffusion  of  the  Increased  pressure  In  the  shock 

front  In  the  downward  direction.  This  in  turn  causes  the  overpressure 
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profile  given  above  to  persist  long  beyond  the  times  given  by  the  t 
decay  rate. 

Additional,  detailed  Information  on  this  effect  may  be  found  In 
(22:4-103,7),  and  (45:231-243). 
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Example  Problem 

The  use  of  the  preceedlng  equations  to  predict  the  conditions  In  the 
path  of  an  aircraft  In  flight  Is  shown  In  the  following  example.  The 
overpressure,  time  of  arrival,  dynamic  pressure  and  gust  loading  are 
calculated  for  an  aircraft  flying  horizontally  directly  away  from  the 
burst  at  a  ground  speed  of  400  miles  per  hour  (587  feet/sec.).  When  the 
shock  hits  the  plane.  It  Is  at  an  altitude  of  6,000  feet  and  a  ground 
range  of  20,000  feet.  The  nuclear  burst  Is  one  megaton  and  occurs  at  an 
altitude  of  3000  feet. 

The  time  of  arrival  Is  given  by  equation  V-1; 

Scaled  Slant  Range,  R  *  (Gnd  Range*  +  Altitude*)  /Yield 
«  (20,000*  +  6000*)  /(lOOO)  -  2,088  feet 

Time  of  Arrival  =  (.543-21 .185(2. 088 )+361. 81(2. 088)*+2383(2.088)*)10 

(1  +  2.04797(2.088)  +2.687^.088)^) 

>  13.7  seconds 

In  order  to  find  the  overpressure  at  altitude,  the  height  of  the 
triple  point  (HTP)  and  the  overpressure  on  the  ground  must  be  known. 

These  two  parameters  can  be  found  from  graphs  in  EM-1  or  from  numerical 
relations  such  as  those  given  by  Potackl.  The  values  given  here  are  from 
EM-1. 

Height  of  the  Triple  Point,  HTP  *  8000  feet 

Overpressure  on  the  ground,  dP  *  4  psl. 

Fractional  Distance  up  Mach  stem  >■  x  >  Height  of  Target/HTP 


6000/8000  -  .75 


From  equation  V-5: 


F  =  .993666  +  4.72762(.75)  -  15.5804(.75)  -  21 .66193( . 75) 

+  127.5835(.75)  -  155.5353( .75)  +  59.9366(.75) 

=  .66 


Overpressure  at  altitude  =  (Overpr.  Gnd)(Tgt.  Press)F/14.7 

=  4(. 66)11. 78/14. 7 
*  2.1  psl. 

The  following  parameters  may  be  found  from  U.  S.  Standard  Atmospheres: 
Target  ambient  air  pressure  °  11.78  psl. 

Target  ambient  air  density  =  .001988  slugs/ft^ 

Target  ambient  air  sound  speed  =  1094  feet/sec. 

Airspeed  of  shock  -  _ 5(1094)2.1 _ 

7(11.78)(1+6(2.1)/7(11.78)) 

■  130  feet/sec. 

Effective  Radius  -  (.574  +  .317(3.)  +  .0000573(3. )*)(20, 000) 

“  30,500  feet 

Angle  of  Shock  ■  T  ®  Arctan(Height  of  Target/Eff.  Radius) 

■  arctan( 6000/30500)  “  11.13  degrees 


Component  of  shockspeed  perpendicular  to  gnd  ~  Ua  ~  UslnT 
=  130sln(11.13)  ■  25  feet/second 
Component  of  shockspeed  parallel  to  gnd  ■>  Ub  ~  Ucos  T 
“  130cos(11.13)  *  128  feet/second 


.-A'-V 
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New  Airspeed  of  Airplane  =  (Ua*  +  (Vj^  _  m,)*) 
-  (25*  +  (587-128)*)  =  460  feet/sec. 


New  angle  of  attack  =  a  =  arcsln(Ua/V) 

=  arcsln(25/460)  =3.1  degrees 

Dynamic  Pressure  =  ^  PV* 

=  ^(.002235slugs/ft*)(460)*/1441n*/ft* 
=  1.6  psl. 


Comparison  With  Conventional  Predictions. 

In  order  to  estimate  the  effect  equation  V-7  would  have  as  compared  to 
conventional  techniques,  a  one  psl.  Iso-overpressure  contour  was  plotted 
for  a  500  Klloton  burst  with  a  height  of  burst  of  3175  feet  (SHOB  =  400 
feet).  The  results  for  calculations  based  upon  equations  IV-6  and  IV-7 
and  conventional  techniques  are  plotted  In  figure  21.  It  Is  clear  that 
the  prediction  method  based  upon  equations  IV-6  and  IV-7  Is  substantially 
different  than  the  traditional  methods  of  prediction.  The  overpressure 
at  altitude  In  the  traditional  method  Is  a  function  only  of  the 
overpressure  on  the  ground  and  the  pressure  at  altitude.  The  proposed 
predictions  are  also  a  function  of  the  fractional  distance  up  the  mach 
stem.  There  Is  a  large  area  Just  below  the  triple  point  that  Is  much 
less  dangerous  than  previously  believed.  Likewise,  nearer  the  ground, 
the  danger  area  extends  much  farther  out  than  before.  In  this  case  of 
the  400  foot  scaled  height  of  burst,  the  danger  area  extends  out  nineteen 
percent  farther  than  the  traditional  prediction  at  3000  feet  altitude. 


Yield  =  500  KT 


Ground 


^  S' 


The  Reflect-4  Codes  have  also  predicted  the  lines  of  constant 
overpressure.  Two  plots  are  shown  as  figures  22  and  23.  On  these 
figures,  the  predictions  given  by  equation  IV-7  for  one  psi.  are  shown. 
For  the  one  megaton  burst  at  3000  feet,  the  two  contours  show  good 
agreement  on  the  lower  half  of  the  contour.  The  discrepancy  between  the 
two  predictions  Is  believed  to  be  caused  by  the  codes'  values  for  the 
angle  tff  In  this  region.  The  Reflect  code  predicts  a  higher  angle  for 
in  the  far  regions  than  other  predictions  (See  (53;figure  6)).  This 
higher  value  for  ^  translates  Into  a  higher  overpressure  In  the  Hach 
stem.  This  also  occurs  In  the  second  plot  for  a  one  megaton  burst  at 
5000  feet.  As  predicted  earlier,  the  magnitude  of  the  overpressure  bulge 
predicted  by  equation  V-5  Is  too  great  for  heights  of  burst  greater  than 
300  scaled  feet.  In  this  case,  the  magnitude  of  the  peak  overpressure  Is 
about  15  percent  too  large. 
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VI.  Conclusions  and  Recommendations 


Summary 

This  paper  has  used  current  shock  theories  augmented  by  empirical 
results  and  new  relations.  In  order  to  trace  the  trajectory  of  the  triple 
point  for  spherical  geometries.  Spherical  shock  configurations  were 
assumed  to  be  Instantaneously  modeled  by  planar  shocks.  This  assumption 
allowed  an  analytical  solution  for  strong  shocks.  Curve  fits  to 
laboratory  data  were  used  when  the  Mach  number  of  the  Incident  shock 
dropped  below  2.  New  relations  were  used  for  the  situation  where  the 
triple  point  Is  above  the  burst  height. 

Empirical  relations  were  developed  In  this  paper  for  the  shape,  time 
of  arrival,  and  overpressure  In  the  Mach  stem.  These  relations  used  data 
from  measurements  taken  from  free-alr  nuclear  detonations  as  well  as  from 
high-explosive  charges  of  varying  sizes.  The  relations  developed  here 
allow  the  calculation  of  the  overpressure  at  any  altitude  In  the  mach 
shock  If  the  triple  point  Is  at  least  twice  as  high  as  the  height  of 
burst. 

The  iso-overpressure  contours  given  by  the  relations  developed  here 
are  markedly  different  from  conventional  calculations.  They  are, 
however,  borne  out  by  the  limited  nuclear  data,  and  the  Reflect-4  Code 
calculations.  The  results  Indicate  that  near  the  ground,  the 
Iso-overpressure  contours  extend  much  farther  away  from  the  blast  than 
previously  believed.  Similarly,  the  contours  near  the  triple  point  are 
much  closer  to  the  burst  point  than  before. 


Discussion 


The  finding  of  a  significant  departure  from  traditional  predictions 
has  several  implications  for  the  military  planner  and  targeteer.  The 
appearance  of  an  area  of  Increased  overpressure  means  that  the  vicinity 
of  the  surface  of  the  ground  is  more  deadly  than  present  predictions 
show.  In  the  Reflect  studies  for  the  Defense  Nuclear  Agency,  it  was 
determined  that  a  pilot  in  the  Mach  stem  had  a  choice  between  climbing  to 
an  altitude  of  lower  density  (and  thus  lower  overpressure)  or  diving  with 
the  Intent  of  increasing  his  velocity  away  from  the  blast  (53:81).  The 
findings  here  indicate  that  he  should  not  decrease  his  altitude.  Diving 
will  not  only  raise  the  ambient  air  pressure  (and  thus  the  overpressure) 
but  will  put  the  plane  closer  to  the  increased  overpressure  region  as 
well. 

An  Increase  in  the  overpressure  above  ground  means  that  it  is  possible 
to  optimize  the  height  of  burst  to  hit  targets  at  some  given  distance  and 
altitude.  In  the  Base  Escape  problem,  there  are  many  airplanes  trying  to 
escape  a  likely  nuclear  blast.  The  probable  distribution  of  these  planes 
in  space  and  time  can  be  determined  and  targeted  to  destroy  the  largest 
number.  Similarly,  large  buildings  can  be  targeted  to  place  the  largest 
dynamic  pressure  at  some  height  to  inflict  the  maximum  damage.  It  is 
Important  to  recognize  that  maximizing  the  overpressure  on  the  ground  is 
not  the  same  thing  as  maximizing  the  overpressure  on  some  structure  of  a 
given  height.  The  'knee'  in  the  overpressure  curves  corresponds  to  the 
maximum  distance  that  regular  reflection  can  be  projected.  At  this 
distance  the  mach  stem  is  not  yet  formed  so  that  the  overpressure  against 
the  side  of  the  building  is  Just  the  free  air  overpressure  and  not  the 
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overpressure  on  the  ground. 

Recommendations 

In  order  to  develop  a  more  efficient  and  reliable  model  of  the 
conditions  to  be  found  in  the  Mach  stem,  two  tasks  need  to  be 
accomplished.  First,  there  must  be  more  measurements  made  In  the  Mach 
stem  region  above  the  surface.  This  will  give  a  better  estimate  of  the 
magnitude  of  the  overpressure  Increase  for  varying  heights  of  burst  and 
the  way  It  propagates  In  space.  This  data  is  crucial  to  test  any 
formulation  purporting  to  predict  conditions  In  the  Mach  stem.  At 
present,  the  data  does  not  cover  a  wide  range  of  burst  heights  and 
distances  well  or  consistently. 

Secondly,  some  type  of  analytical  technique  must  be  developed  to 
predict  the  overpressure  In  the  Mach  stem  at  distances  far  below  the 
triple  point.  At  present,  theory  predicts  only  the  overpressure  at  some 
Infinitesimal  distance  below  the  triple  point.  Knowing  the  overpressure 
variation  with  altitude,  one  must  be  able  to  predict  the  manner  in  which 
the  shock  will  propagate  In  space.  The  shock  will  clearly  deform.  This 
is  Important  to  the  calculation  of  gust  loading  of  airfoils.  The 
overpressure  bulge  will  tend  to  diffuse  and  equalize  the  pressure  in  the 
shock  throughout  the  front.  This  means  the  magnitude  of  the  increase  in 
overpressure  will  change  with  time.  The  only  way  known  to  this  author  of 
predicting  the  manner  in  which  this  kind  of  system  will  propagate  Is  the 
Whltham  shock  ray  theory.  (67:235).  This  theory  has  been  used  with  some 
success  In  predicting  the  effects  of  broken  terrain  on  shock  waves 
(22:5-100).  The  result  of  this  would  be  a  more  comprehensive  compilation 
of  curves  by  which  one  could  estimate  the  magnitude  and  location  of  the 


overpressure  bulge  at  any  time  and  for  any  height  of  burst.  This 
information  would  be  of  great  value  to  military  planners  who  must  plan 
and  protect  against  such  possibilities. 

Present  methods  of  predicting  Mach  stem  parameters  are  clearly 
Inadequate.  Until  more  data  is  available  and  properly  analysed,  Mach 
stem  models  will  remain  Inaccurate,  unreliable,  and  limited  in  usefulness 
at  best. 
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APPENDIX  A 


Oblique  Shock  Relations 

The  change  In  state  across  a  shock  Is  found  from  consideration  of 
conservation  of  mass,  momentum,  and  energy  across  a  control  volume 
Including  the  shock  front.  This  will  give  three  equations  linking 
temperature,  pressure,  density,  and  velocity  across  the  front.  The 
equation  of  state  for  air  Is  a  fourth  equation  to  link  parameters 
together  on  the  same  side  of  the  shock.  These  equations  have  been 
combined  In  many  forms  to  obtain  one  set  of  parameters  In  terms  of  any  of 
the  others.  Several  of  the  most  useful  In  the  analysis  of  this  problem 
have  been:  (56:160) 


Jli  “  UL  H  »sttv»0  -  XzL 
^  y+i  ^  y+i 

^  ^  -  V,8ln0  =  tan(  <f> ) 

Pi  Vn2  V^slirr^-  d )  tanC0'-  6  ) 

Vti  =  Vt2 

Vicos0  =  V2Cos(^-i  ) 

Where:  Vi,  are  the  velocities  of  the  approaching  and  exiting  air 
Vti,  Vt^  are  the  velocity  components  tangent  to  the  shock 
VOi,  Vn2  are  the  velocity  components  normal  to  the  shock 
Pl,  P2  are  the  densities  on  each  side  of  the  shock 
Is  the  angle  of  Incidence 
6  Is  the  angle  of  deflection  of  the  Incident  air 
Pl,  P2  are  the  pressures  on  each  side  of  the  shock 


M  Is  the  mach  number  of  the  shock 


y  is  the  ratio  of  specific  heats,  for  air  equal  to  1.4 


Also  of  use  are  the  Ranklne-Hugoniot  equations.  This  set  set  of 
equations  as  commonly  written  are  derived  from  the  above  relations  for  a 
normal  shock.  They  express  the  conditions  across  the  shock  front  in 
terms  of  the  similar  property  ahead  of  the  shock  and  the  overpressure. 
Overpressure  is  the  easiest  property  of  a  shock  wave  to  measure  and 
probably  the  best  predicted  property  of  an  incident  shock  from  a 
detonation. 


dP  is  the  pressure  difference  across  the  shock 
All  other  symbols  are  the  same 


APPENDIX  B 

Mach  Stem  Overpressure  Data 


Nuclear  Data 


Cannlster  Amb.  Press 
Psl. 

Gnd  Dlst 
Feet 

Vert.  Alt 
Feet 

Overpress 

Psl. 

8 

Teapot  4  (30: 

11.60 

27) 

6069’ 

6769' 

3.78 

3.77 

9 

10.98 

5386’ 

8319’ 

3.00 

3.914 

10 

10.37 

5172’ 

9819’ 

1.90 

4.469 

6 

Ivy  King  (31:27) 

7.17  69081’ 

19550’ 

.27 

58.62 

7 

11.46 

39946’ 

7050’ 

.875 

29.93 

9 

10.30 

40109’ 

10000’ 

.805 

30.48 

10 

10.42 

39150’ 

9700’ 

.805 

30.36 

11 

9.80 

40570’ 

11400’ 

.72 

31.33 

7 

Upshot-Knothole 

11.70 

9  (34:25) 

34568’ 

6550’ 

.255 

29.98 

10 

9.46 

52321’ 

11925’ 

.11 

47.39 

15 

11.80 

27764’ 

6325' 

.37 

23.65 

18 

11.08 

34327’ 

7875' 

.235 

29.76 

19 

10.84 

36380’ 

8400’ 

.22 

31.83 

20 

10.87 

38734’ 

8350’ 

.205 

33.93 

1 

Snapper  5  (33 

5.89 

:25,20) 

26157’ 

32900’ 

.141 

26.99 

2 

6.14 

16754’ 

25150’ 

.175 

20.52 

3 

4.83 

19579’ 

32000’ 

.128 

25.81 

4 

7.66 

7235’ 

14950' 

.425 

11.58 

nnlster  Amb.  Press 

5  4.50 

Gnd  Dlst 
8676* 

Vert.  Alt 
27750’ 

Overpress 

.147 

.  Ta 

23.37 

6 

6.43 

9334’ 

18500’ 

.310 

16.09 

8 

5.23 

17096' 

25250' 

.168 

24.34 

10 

8.50 

7559* 

12760’ 

.604 

9.59 

12 

11.01 

3725’ 

5500’ 

2.69 

3.05 

13 

8.94 

1345’ 

9100’ 

.817 

6.28 

15 

8.60 

10515’ 

13350’ 

.486 

11.68 

16 

8.34 

19727’ 

19050’ 

.307 

19.17 

1 

Snapper  8  (33: 

6.45 

26,27,33) 

11135’ 

21950’ 

.276 

17.13 

4 

8.94 

9053’ 

13760’ 

.624 

10.17 

6 

7.52 

14549’ 

18150’ 

.314 

16.55 

7 

5.54 

14717’ 

25630’ 

.207 

22.13 

8 

6.63 

25815’ 

21300’ 

.205 

26.59 

9 

9.22 

17596’ 

12960’ 

.370 

15.73 

10 

12.19 

8016’ 

5240’ 

1.38 

5.67 

11 

11.79 

9040’ 

6210’ 

1.033 

6.73 

14 

9.24 

4038’ 

12900’ 

.800 

6.95 

15 

11.38 

3662’ 

7200’ 

3.235 

2.74 

16 

8.85 

11479’ 

14000’ 

.488 

11.78 

A 

Plumbob-Owens 

11.00 

(25:23) 

9605’ 

2900’ 

.75 

2.42 

B 

10.85 

9281’ 

3800’ 

.94 

2.75 

D 

9.80 

7200’ 

7000’ 

.65 

2.83 

Redwing  (32) 

Data  Available  In  Reference  32 


Notes:  All  Owens  Alt.  are  height  above  ground 
All  ground  distances  calculated. 


TNT  Bursts  Above  Ground 


HOB 

Gnd  Dlst. 

HTP 

Height  of  Gage 

Overpressure 

2.31' 

13.67' 

3.4' 

0.0' 

1.43-1.47 

•• 

■1 

•• 

.91’ 

2.18 

It 

II 

•I 

1.9' 

2.18 

II 

II 

•1 

2.8’ 

1.8 

II 

19.67' 

7.0' 

0.0' 

1.29-1.31 

■1 

•I 

M 

.91’ 

1.81 

II 

II 

•  1 

2.8' 

1.54 

■I 

•I 

•I 

5.75’ 

1.40 

•I 

ti 

tl 

6.75’ 

1.2 

3.42' 

19.67' 

3.5' 

0.0' 

1.54-1.62 

•1 

II 

II 

.91’ 

2.02 

M 

•1 

II 

2.8’ 

1.64 

11 

22.67' 

4.5' 

0.0’ 

1.43-1.56 

•• 

•1 

II 

.91’ 

1.98 

•1 

•1 

II 

2.8’ 

1.58 

4.33' 

16.67' 

1.4' 

0.0' 

1.84-2.04 

•1 

II 

II 

.91’ 

2.22 

fi 

19.67' 

1.8' 

0.0’ 

1.74-1.92 

•• 

ti 

II 

.91’ 

2.44 

M 

22.67' 

2.5’ 

0.0' 

1.67-1.81 

•I 

M 

.91' 

2.13 

Notes:  All  bursts  are  \  lb.  of  TNT 

All  data  from  Stoner  (62:31-34) 

Free  air  data  In  same  tests  show  a  Std.  Dev.  of  6.5  percent 


Project  Dipole  West  Detonations 


Shot 

HOB 

Gnd  Dlst 

HTP 

Height  of  Gage 

Overpresssure 

7 

25.46’ 

60’ 

7’ 

0’ 

27.5 

7 

25.46’ 

60’ 

7’ 

3’ 

36.0 

8 

24.45’ 

60’ 

7’ 

0’ 

31.5 

8 

24.45’ 

60’ 

7’ 

3’ 

34.8 

9 

15.15’ 

40’ 

5* 

0’ 

79 

9 

15.15 

40’ 

5’ 

3’ 

80 

9 

15.15’ 

60' 

15+’ 

0’ 

35.6 

9 

15.15’ 

60’ 

15+’ 

3’ 

33.5 

9 

15.15’ 

60’ 

15+’ 

10’ 

52.4 

10 

14.92’ 

40’ 

7’ 

0’ 

67.5 

10 

14.92’ 

40’ 

7’ 

3’ 

68 

10 

14.92’ 

60’ 

15’ 

0’ 

30.8 

10 

14.92’ 

60’ 

15’ 

3’ 

30.0 

10 

14.92’ 

60’ 

15’ 

10’ 

43 

Shot  7  data  (37:143) 

Shot  8  data  (37:151) 

Shot  9  data  (37:158-9) 

Shot  10  data  (37:165-6) 

All  shots  were  1080  lbs  of  penollte 
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Planar  Shock  Data 


Angle 

from  Wedge  Fract 

Overpr  Gnd  Overpr  at  Gage 

Fract 

Case 

10  Mach 

No.  1.66 

Wedge  Angle 

=40 

=3.3 

1 

.303 

6.33 

7.48 

1.2 

Case 

11  Mach 

No.  1.90 

Wedge  Angle 

-40 

=3.8 

1 

.263 

8.92 

12.04 

1.34 

3.5 

.92 

9.35 

8.41 

.89 

Case 

12  Mach 

No.  2.41 

Wedge  Angle 

=40 

=4.4 

1 

.227 

12.86 

18.67 

1.45 

3.5 

.795 

13.44 

11.62 

.76 

Case 

13  Mach 

No.  2.87 

Wedge  Angle 

-40 

-4.5 

1 

.222 

24.20 

29.86 

1.24 

3.5 

.77 

24.58 

18.78 

.76 

Notes:  Data  from  ref.  17  table  3,  figures  9  and  10. 

Ground  overpressures  are  Interpolated  from  measurements 
1.33  cm  apart 


APPENDIX  C 

Origin  of  the  Triple  Point  Data 


HIGH  EXPLOSIVE  TESTS 


HOB 

Shock  Strength 

Gnd  Range 

Extreme  Angle 

.5' 

.0085 

.417' 

39.80 

1.0' 

.02280 

.833' 

39.80 

1.25' 

.03306 

1.000' 

39.667 

1.50' 

.04858 

1.250' 

39.80 

2.00' 

.07842 

1.667' 

39.80 

3.00' 

.17464 

2.500' 

39.80 

4.00' 

.30225 

3.458' 

40.85 

5.00' 

.43867 

4.458' 

41.72 

6.00 

.56994 

5.916' 

44.60 

■2a. 

i 


Data  From  Bryant  (12:17) 

The  above  data  was  collected  from  I  Ib.  TNT  bursts  over  hardpacked 
earth  and  dry  sand.  It  should  be  noted  that  while  the  trajectories  of 
the  triple  points  were  different  for  earth  and  sand,  the  origins  of  the 
triple  points  are  the  same  within  the  limits  of  measurement  (12:16). 
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NUCLEAR  EXPLOSIVE  TESTS 


Shot 

HOB 

Gnd  Range 

OverPress. 

Strength 

Extreme  Angle 

Teapot  4 

500' 

550' 

200 

.07 

47.8 

Teapot  6 

500' 

450' 

162 

.08 

42.0 

Teapot  8 

500' 

340' 

450 

.03 

34.2 

Teapot  12 

400' 

325' 

450 

.03 

39.1 

Tumbler  1 

793' 

497' 

2.4 

.86 

32.07 

Tumbler  2 

1109' 

749' 

4.97 

.747 

34.03 

Tumbler  3 

3447' 

2030' 

5.77 

.71 

30.63 

Priscilla 

700' 

550' 

366 

.038 

38.2 

Tumbler  Data  (2:56) 
Priscilla  Data  (63:62) 


The  ground  range  to  the  OTP  in  some  cases  is  estimated.  It  should  be 


noted  that  different  measurement  techniques  often  give  conflicting 
results.  Because  of  the  precursor,  gage  data  is  unreliable  and 
photographs  may  be  obscured  by  dust  and  debris. 


APPENDIX  D 


Near  Ground  Overnressure  Data 


Terrain  Gnd  Dlst 

Height 

Overpr 

Terrain 

Gnd  Dlst 

Height 

Overpr 

Teapot  12 
Water  2000' 

(55:46- 

0' 

■51) 

17.4 

Water 

2500' 

O' 

11.8 

2000’ 

3’ 

20.1 

2500’ 

3’ 

13.2 

2000' 

10' 

18.1 

2500’ 

10' 

12.9 

Water  3000' 

O' 

8.76 

2500' 

25' 

13.7 

3000' 

3’ 

10.5 

2500' 

40’ 

11.2 

Desert  2000' 

O' 

16.9 

Asphalt 

2500' 

0' 

6.6 

2000' 

3' 

18.6 

2500' 

3' 

8.5 

2000' 

10' 

21.9 

2500' 

10’ 

6.32 

2500' 

25' 

6.9 

2500' 

40' 

6.8 

Tumbler  4 
7300' 

(1:44) 

O' 

2.17 

5000' 

O' 

3.82 

7300' 

10' 

2.17 

5000' 

10' 

3.85 

7300' 

50' 

2.17 

5000' 

50' 

3.88 

Upshot  Knothole  9 
Classified 

(45:64, 

65,66) 

APPENDK  E 


Nuclear  Explosion  Parameters 

The  data  presented  here  Is  a  summary  of  the  nuclear  explosions  that 
were  of  use  in  this  analysis.  It  must  be  recognized  that  all  of  the  data 
presented  here  Is  not  unquestioned,  particularly  the  given  yields. 
Different  references  may  report  substantially  different  yields,  heights 
of  burst,  and  surface  pressures. 


Title 

Yield 

HOB 

Seal  HOB 

Press.  Gnd 

Alt.  Gnd 

Ref. 

U-K  9 

26 

2423' 

763' 

13.05 

4191' 

(10:78) 

Tumbler  1 

1.05 

793’ 

747' 

13.26 

— 

(2:202) 

Tumbler  2 

1.15 

1109' 

995* 

12.73 

— 

(2:202) 

Tumbler  3 

30 

3447' 

1012’ 

11.17 

— 

(2:202) 

Tumbler  4 

19.6 

1040’ 

363' 

12.72 

— 

(2:202) 

Priscilla 

36.6 

700' 

201.7' 

13.12 

3078' 

(13:29) 

Owens 

9.7 

500' 

222' 

12.60 

— 

(13:29) 

Teapot  4 

43 

500' 

135' 

12.60 

4491' 

(46:17) 

Teapot  6 

7.76 

500' 

240’ 

12.71 

4  245' 

(46:17) 

Teapot  8 

14.23 

500' 

195' 

12.58 

4309' 

(46:17) 

Teapot  12 

22.0 

400' 

137’ 

13.18 

3077' 

(46:17) 

Ivy  King 

550 

1500' 

183’ 

14.68 

0.0' 

(31:20,45) 

Snapper  5 

11.7 

300' 

126' 

12.78 

4200' 

(33:35) 

Snapper  8 

14 

300’ 

118’ 

12.68 

4200’ 

(33:35) 

APPENDK  F 
Program  Description 

A  program  was  written  using  the  techniques  described  in  this  paper  to 
predict  the  conditions  in  the  path  of  the  mach  stem.  It  is  written  in 
BASIC  for  a  Texas  Instruments  computer.  The  language  is  a  standard  BASIC 
except  for  the  way  that  subroutines  are  called.  Subroutines  are  CALLed 
just  as  in  FORTRAN  using  a  parameter  list  to  pass  variables.  Exits  from 
a  subroutine  are  accomplished  using  a  SUBEXIT  or  a  SUBEND  statement. 

This  program  is  intended  to  be  a  subroutine  of  a  much  larger  program  for 
general  predictions  from  nuclear  bursts.  The  program  is  called,  passing 
the  values  for  the  scaled  ground  range  (SGR),  scaled  height  of  burst 
(SHOB),  scaled  height  of  the  target  (SHT),  the  ambient  pressure  at  the 
ground  (PG) ,  and  the  pressure  at  the  target  (PT) .  The  parameters 
returned  are  the  dynamic  pressure  (Q),  the  overpressure  at  the  target 
(OPT),  and  the  angle  of  the  dynamic  impulse  (TAU). 

There  are  four  main  subroutines  to  the  program.  They  are  TPLPT, 
GNDOVERPR,  TGTOVERPR,  and  DYNAMICPR.  Of  these,  TPLPT,  which  calculates 
the  triple  point  trajectory,  is  the  largest.  TPLPT  has  a  main  calling 
routine  that  starts  at  a  ground  range  in  the  regular  reflection  region 
and  works  its  way  outward,  calculating  the  value  of  ^  and  the  height  of 
the  triple  point  as  it  goes.  The  ground  range  is  stepped  in  Increments 
of  five  percent  of  the  scaled  height  of  burst  until  the  ground  range  of 
the  target  is  exceeded.  On  this  last  iteration,  the  height  of  the  triple 
point  at  the  target  ground  range  is  found  by  interpolation  between  the 
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last  two  values. 


The  subroutine  Is  based  upon  knowing  the  height  of  the  triple  point  at 
some  ground  range.  From  this  Information,  the  slant  range  to  that  point 
can  be  found.  Subroutine  OVERPR  Is  called  to  find  the  overpressure  and 
mach  number  of  the  Incident  shock  at  that  point.  The  geometry  of  the 
problem  allows  determination  of  the  angle  of  Incidence.  With  this 
Information  the  subroutine  PSI  Is  called  to  find  the  angle  at  which  the 
triple  point  rises.  Knowing  the  step  size  and  the  angle,  the  location  in 
space  of  the  next  point  can  be  found.  Once  two  Iterations  are  done,  the 
triple  point  Is  found  at  the  next  point  by  Integration  of  equation  IV-1. 
The  subroutine  INTEG  Is  used  to  integrate  the  equation  using  Simpson's 
rule . 

Subroutine  PSI  uses  the  techniques  outlined  In  Section  IV  of  this 
paper.  It  Is  divided  Into  three  parts.  One  calculates  the  angle  0  «if 
the  triple  point  Is  above  the  height  of  burst.  The  other  two  are  used 
only  when  the  triple  point  Is  below  the  height  of  burst  and  are 
differentiated  by  whether  the  mach  number  of  the  Incident  shock  Is 
greater  or  less  than  2.5. 

The  overpressure  on  the  ground  is  found  by  use  of  the  subroutine 
GNDOVERPR.  This  routine  uses  the  Horizons  Technology  version  of  the  EM-1 
curves.  Input  to  the  program  consists  of  the  scaled  ground  range  at 
which  the  overpressure  Is  desired  and  the  scaled  height  of  burst.  The 
overpressure  on  the  ground  Is  then  returned  after  being  scaled  to  the 
ground  ambient  pressure.  This  overpressure  Is  then  Input  Into  the 
subroutine  TGTOVERPR  to  find  the  overpressure  at  altitude.  This  Is  found 
by  scaling  to  the  pressure  of  the  target  and  multiplying  by  a  factor. 


This  factor  is  described  In  the  text  and  describes  the  pressure  variation 
with  height. 

The  last  of  the  routines  Is  the  calculation  of  the  dynamic  pressure. 
This  is  done  by  the  use  of  the  Rankine-Hugoniot  equations.  The  angle  of 
the  impulse  Is  also  found.  This  Is  found  through  the  estimate  of  the 
mach  stem  radius  as  described  In  the  text.  This  angle  Is  used  to  find 
the  effective  angle  of  attack  In  the  gust  loading  of  an  airfoil. 


100  ! 

110  !***  SUBROUTINE  MACH  *** 

120  !***  CALCULATES  THE  CONDITION  INSIDE  *** 

130  1***(;THE  MACH  STEM  IN  A  NUCLEAR  BRURST  *** 

140  !  *****'if***********it****^‘'tt******‘fi^i*******it* 

150  SUB  MACH(SGR,SHOB,PT,PG,SHT,Q»OPT,TAU) 

160  !  SGR - SCALED  GROUND  RANGE  (FEET) 

170  !  SHOE— SCALED  HEIGHT  OF  BURST  (FEET) 

180  !  SHT - SCALED  HEIGHT  OF  TARGET  (FEET) 

190  !  HTP - HEIGHT  OF  THE  TRIPLE  POINT  (FEET) 

200  I  OPG - OVERPRESSURE  ON  THE  GROUND  (PSI) 

210  I  OPT - OVERPRESSURE  ON  THE  TARGET  (PSI) 

220  !  PG - ATMOS.  PRESSURE  ON  THE  GROUND  (PSI) 

230  !  PT - ATMOS.  PRESSURE  AT  THE  TARGET  (PSI) 

240  !  Q - DYNAMIC  PRESSURE  (PSI) 

250  CALL  TPLPT(9GR,SH0B,HTP) 

260  IF  HT§HTF  THEN  290 

270  PRINT  "TARGET  ABOVE  THE  MACH  STEM" 

280  SUBEXIT 

290  CALL  GND0VERPR(SH0B,S6R,0PG) 

300  CALL  TGTOVERPR(OPG,OPT,PG,PT,HT,HTP) 

310  CALL  DYNAMICPR(Q,OPT,PT,TAU,96R,SROB) 

320  SUBEND 

330  ! ***************************************** 

340  !***  *** 

350  I***  TRIPLE  POINT  ROUTINE  *** 

360  I***  *** 

370  ! ***************************************** 

380  SUB  TPLPT(SGR,SHOB,YT) 

390  I  X - GROUND  RANGE/HOB 

400  !  TH - 90-ANGLE  OF  INCIDENCE 

410  !  FSl - ANGLE  AT  WHICH  THE  TRIPLE  POINT  RISES 

420  t  M - MACH  NUMBER  OF  THE  INCIDENT  SHOCK 

430  !  DP - OVERPRESSURE  OF  INCIDENT  SHOCK 

440  !  YT - TRIPLE  POINT  HEIGHT/HOB 

4  50  X-0.75 
460  Xl-SGR/SHOB 
470  GAMMA-1.4 
480  HOB-SHOB/3.28 

490  I ***************************************** 

500  I***  TRIPLE  POINT  *** 

510  !***  MAIN  PROGRAM  ROUTINE  *** 

520  I***  *** 

530  I ***************************************** 

540  CALL  OVERPR(HOB,X,YT,M,DF) 

550  TH-ATN((1-YT)/X) 

560  CALL  PSI(H0B^,YT,TH,PSI1,GAMMA,M,DP) 

570  X««+.05 


580  YT1-YT+.05*TAN(PSI1) 

590  TH-ATN((1-YT1)/X) 

600  CALL  0VERPR(H0B^,YT1,M,DP) 

610  CALL  PSI(H0B^,YT1,TH,PSI2,GAMMA,M,DP) 

620  CALL  INTEG(PSI1,PSI2,YT,X) 

630  IF  X§!(1  THEN  670 

640  YT-YT-10*(X-X1)/(YT-YT0) 

650  YT-YT*SH0B 
660  SUBEND 
670  X^+.05 
680  YTO-YT 
690  GOTO  540 

710  !***  *** 

720  !***  SUBROUTINE  INTEGRATE  *** 

730  !***  *** 

740  ! ***************************************** 

750  SUB  INTEG(PSI1,PSI2,YT^) 

760  1  SUBROUTINE  USES  SIMPSONS  RULE  TO  INTEGRATE  AND  PINT  YT 
770  YT-YT+.05*(TAN(PSI0)+4*TAN(PSI1)+TAN(PSI2))/3 
780  PSI0-PSI2 
790  SUBEND 

800  I  ***************************************** 

810  !***  *** 

820  !***  SUBROUTINE  OVERPR  *** 

830  !***  *** 

840  ] ***************************************** 

850  SUB  OVERPR(HOB,X,YT,M,DP) 

860  I  SUBROUTINE  USES  HORIZONS  TECHNOLOGY  ROUTINE  TO  FIND  FREE  AIR  OVERPRESSURE 

870  1  M  IS  FOUND  USINC  RANKINE-HUGONIOT  EQUATIONS 

880  R1-H0B*(X*X+(1-YT)*(1-YT))C.5 

890  DP-EXP( . 19*(LOG(R1/1000) )c2-l .5*LOG(Rl/1000)-. 1 ) 

900  M-(1+(6*DP)/(7*14.7))C.5 
910  SUBEND 

920  I ***************************************** 

930  1***  *** 

940  I***  SUBROUTINE  PSI  *** 

950  1***  *** 

960  i ***************************************** 

970  I***  *** 

980  !***  MACH  NUMBERS  V  2.5  *** 

990  I***  *** 

1000  ! **************************************** 

1010  SUB  PSI(HOB^,YT,TH,PH,GAMMA»M,DP} 

1020  IF  TH§0  THEN  1910 
1030  IF  M§2.5  THEN  GOTO  1450 
1040  PHL-.8*PH 
1050  PH-PI/180fPHL 


1060  PHH-PI/2 

1070  K1-(GAMMA+1)/(GAMMA-1) 

1080  K2-2*GAMMA/(GAMMA+1) 

1090  TT=TAN(TH) 

1100  TP-TAN(PH) 

1110  K1@-(K1*M*M*COS(TH+PH)c2)/(2/(GAMMA-1)+M*M*COS(T1H-PH)c2) 

1120  K2@-K2-l/(Kl*M*M*C0S(TftfPH)c2) 

1130  CT-1/TAN(TH) 

1140  CP-1/TAN(PH) 

1150  C0T-1/TAN(P1H-TH) 

1160  E-((K2*M*M*COS(PH)-1/K1)/(K2*M*M*COS(T1H-PH)-1/K1))c2 

1170  A1-(CP/K1@-(TT+C0T/K1@)/(1-TT*C0T/K1@))/(1+CP*((TT+C0T/K1@)/(1-TT*C0T/K1@)) 
/Kl@) 

1180  A2-(Kl*K2*(l+Kl@*Kie*(TAN(TttfPH))c2))/(Kl@*K2@*(l+Kl*E))-l 
1190  A3»(E-l)/(GAMMA*(l+Kl@*Kl@*(TAN(PiH-TH))c2)/(Kie*K2@)-E+l) 

1200  DELTA=A1*A1-A2*A3*A3 
1210  IF  (PHH-PHL)SlE-5  THEN  1300 
1220  IF  ABS(DELTA)§. 00001  THEN  1300 
1230  IF  DELTASO  then  1270 
1240  PHL-PH 

1250  PH-PHfMIN(PI/180,(PHH-PH)/2) 

1260  GOTO  1090 
1270  PHH-PH 
1280  PH-(P8+PHL)/2 
1290  GOTO  1090 
1300  SUBEXIT 

1310  f  ***************************************** 

1320  !***  *** 

1330  !***  MACH  NUMBERS  S  2.5  *** 

1340  {***  *** 

1350  ! ***************************************** 

1360  DATA  -679. 4273854,44. 29925276,-.787886744,-. 0075477956, .0004252447,-5.11668 
E-6 , 2 . 06509E-8 ,-7.01 72E-2 1 

1370  DATA  -87 . 33954,-4 .195549 , .8139311 ,-.035075 , .000691789,-6 . 58892E-6 ,2 . 46545E- 
8,-4.73696E-21 

1380  DATA  103. 1200447,-6. 369145,. 026000227,. 00615099,-. 00017623, 1.94109E-6,-7. 76 
64E“9  2  08813E*22 

1390  DATA  18. 05405, 1.215717896,-. 2132145477,. 0088911,-. 0001684565, 1.55164E-6, -5. 
61549E-9,5.64009E-21 

1400  DATA  -651 .37515,58. 00225446 ,-2.07974387 , . 0376534976 ,- . 00035084 , 1 . 47679E-6 ,- 
1  54965E“9  ~6«23181E*21 

1410  DATA  -4. 38396,-. 1800461545,. 018402893,-. 0006043663,. 0000104168,-8. 80233E-8, 
3 . 33744E-10 ,-l . 76475E-22 

1420  DATA  3382. 683563,-269. 85455,8. 74770105,-. 1465537887, .0013232252, -5. 97838E-6 
,1.01429E-8,-4.32411E-22 


1430  DATA  -227 . 2637 , 13 . 76562055 ,-.3163955817 , .0031118296,-5 . 22245E-6 ,-l . 29688E-7 

, 7 . 76104E-10 ,-l . 4341E-21 

1440  DATA  0,0, 0,0, 0,0, 0,0 

1450  T-(PI/2-TH)*180/PI 

1460  E-l/((M*M-l)*7/6+l) 

1470  IF  Eir.l5  THEN  1510 

1480  RESTORE  1360 

1490  DELTA-1  ::  EO-0 

1500  GOTO  1810 

1510  IF  Eir.3  THEN  1550 

1520  RESTORE  1360 

1530  DELTA-. 15  ;;  E0-.15 

1540  GOTO  1810 

1550  IF  Eir.4  THEN  1590 

1560  RESTORE  1370 

1570  DELTA-. 10  E0-.30 

1580  GOTO  1810 

1590  IF  Eir.5  THEN  1630 

1600  RESTORE  1380 

1610  DELTA-. 10  ::  E0-.40 

1620  GOTO  1810 

1630  IF  ET.6  THEN  1670 

1640  RESTORE  1390 

1650  DELTA-. 10  ::  E0-.50 

1660  GOTO  1810 

1670  IF  ET.7  THEN  1710 

1680  RESTORE  1400 

1690  DELTA-. 10  ::  E0-.60 

1700  GOTO  1810 

1710  IF  ET.8  THEN  1750 

1720  RESTORE  1410 

1730  DELTA-. 10  ;;  E0-.70 

1740  GOTO  1810 

1750  IF  EV.9  THEN  1790 

1760  RESTORE  1420 

1770  DELTA-. 10  E0-.80 

1780  GOTO  1810 

1790  RESTORE  1430 

1800  DELTA-1E65  ::  E0-.9 

1810  READ  A(1),A(2),A(3),A(4),A(5),A(6),A(7),A(8) 

1820  READ  B(1),B(2),B(3),B(4),B(5),B(6),B(7),B(8) 

1830  Al-A( 1 )+A<  2 ) *T+A( 3 ) *Tc2+A( 4 )*Tc^+A( 5 ) *Tc4+A( 6 )*Tc5+A( 7 ) *Tc6+A( 8 ) *TC  7 
1840  B1-B(1)+B(2)*T+B(3)*T«2+B(4)*Tc3+B(5)*Tc4+B(6)*Tc5+B(7)*Tc6+B(8)*Tc7 
1850  PHI-A1+(B1-A1)*(E-E0)/DELTA 
1860  PHI-MAX(0,PHI) 

1870  PH-PHI*PI/180 
1880  SUBEXIT 


1890  i ***************************************** 

1900  !***  *** 

1910  1***  NEGATIVE  WEDGE  ANSLES  *** 

1920  !***  *** 

1930  I  ***************************************** 

1940  Yl— YT 

1950  CALL  0VERPR(H0B^,Y1,M1,D2) 

1960  F-((Dm4.7)*(7*14.7+DP)/(14.7*(7*14.7+6*DP)))«;.5 

1970  M1»M1*F 

1980  T-ATN((YT+1)/X) 

1990  PH-2*ATN((SIN(Pl/2+TH)*Ml-SIN(Pl/2-T)*M)/(C0S(Pl/2-T)*M-C0S(Pl/2+TH)*Ml)) 
2000  SUBEND 

2010  ! ***************************************** 

2020  1***  *** 

2030  !***  GROUND  OVERPRESSURE  *** 

2040  !***  *** 

2050  ! ***************************************** 

2060  SUB  GNDOVERPR(SHOB,SGR,OPG) 

2070  1  SUBROUTINE  USES  HORIZONS  TECHNOLOGY  ROUTINE  FOR  GROUND  OVERPRESSURE 

2080  HOB-SHOB/3.28 

2090  GRR-SGR/3.28 

2100  DELTA-ATN(HOB/GRR) 

2110  SLS=SQR(HOB<;2+GRRc2) 

2120  X-LOG(SLS) 

2130  P90».01*EXP(40.3*(SLSc-.295)) 

2140  PTH-.001*EXP(31.3*<SLS«-.2136)) 

2150  DELTAPA-P90-(P90“PTH)*(C0S(DELTA))c2  ( 

2160  IF  SLSV-lOO  THEN  2190 
2170  DELTAPR-DELTAPA 
2180  SUBEXIT 

2190  ALPHA-EXP( .3549*Xc3-6.7133*X<i2+41.468*X-82.819) 

2200  BETA-EXP( .25192*X«;4-5.8741*X<?Sf50.298*Xc2-185.95*X+248.8) 

2210  GAMMA-EXP(.1826*Xc4-4.36786*X<!3+38.6017*Xc2-149.59*X+216.26) 

2220  DELPTAPB-( (C0S(DELTA)« (2*BETA) )*( (SIN(DELTA) )CALPHA)*EXP(GAMMA) 

2230  DELTAPR-DELTAPA+DELTAPB 
2240  DELTAPX-DELTAPR*PG/14.7 
2250  SUBEND 

2260  i ***************************************** 

2270  !***  TARGET  OVERPRESSURE  *** 

2280  ] ***************************************** 

2290  SUB  TGTOVERPR(OPG,OPT,PG,PT,HT,HTP) 

2300  !***  *** 

2310  I***  *** 

2320  X-HT/HTP 

2330  F- .  993766+4 . 72762*X-15 . 580*X«2-21 . 6619*Xc3+127 . 5835*Xc4-155 . 5353*Xc5+59 . 936 
6*Xc6 

2340  OPT-OPG*F*PT/PG 
2350  SUBEND 


2360  ! ***************************************** 
2370  !***  *** 

2380  !***  DYNAMIC  PRESSURE  CALCULATION  *** 
2390  I***  *** 

2400  ! ***************************************** 
2410  SUB  DYNAMICPR(Q,0PT,PT,TAU,9GR,SHT) 

2420  !  Q  IS  FOUND  FROM  RANKINE-HUG ONIOT  EQUATIONS 
2430  Q-2.5*OPT/(7*PT+OPT) 

2440  H«SHOB/100 

2450  RE=(.574+.317*lH-.0000573*Hc2)*SGR 
2460  TAU-ATN(SHT/RE) 

2470  SUBEND 
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